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Computed Potential Energy Surfaces for Chemical Reactions 
Investigator: Dr. Stephen Walch 


In the previous 6 month progress report, calculations of potential energy sur- 
faces (PES’s) for CH3 + OH [1] , CH3O [2] , CH( 2 n) + N2 [3] , and for NH2 + 

NO [4] were described. This work has since been accepted for publication. In 
addition, calculations of the PES for NH2 + 0 have been published [5], and cal- 
culations of the PES for CH 3 + 0 2 [6] have been submitted for publication. ( 
Reprints or preprints of refs. 1-6 are included in the appendix.) 

During the time period from 1 January 1993 to 31 August 1993 several new 
projects have been undertaken. In collaboration with A. Kuppermann ( Cal Tech) 
a new global potential energy surface is being generated for 

0{ z P) + H 2 ^0H + H (1) 

This surface is being fit using the rotated Morse oscillator method, which was 
used to fit the previous POL-CI surface (Walch, Dunning, and Wagner 1980). 

The new surface is expected to be more accurate and also includes a much more 
complete sampling of bent geometries. To date, 350 points have been generated 
for the 3 A' and 3 A” components of the 3 II surface. These 350 points should be 
sufficient to generate the initial fit. The resulting PES will be used in benchmark 
coupled- channel calculations for this system. 

A new study has been undertaken of the reaction 

N + 0 2 - NO + O (2) 

which had been studied previously at Ames (Walch and Jaffe 1987). The new 
studies have focused on the region of the surface near a possible minimum cor- 
responding to the peroxy form of NOO. The previous work of Walch and Jaffe 
had indicated a shallow minimum on the PES for NOO at the CASSCF level, but 
subsequent externally contracted Cl calculations, along a somewhat constrained 
minimum energy path, showed no minimum. More recently, Meredith and Schae- 
fer found a minimum for NOO using CISD gradient methods, in apparent contra- 
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diction to the earlier work of Walch and Jaffe. The new studies use the internally 
contracated Cl methodology of Werner and Knowles ( MOLPRO ), which permits 
much more extensive Cl calculations than were possible in 1987. The new calcu- 
lations show a very shallow minimum, corresponding to N00; however, there is a 
saddle point leading to NO + 0 in close proximity to the minimum, and the bar- 
rier to formation of NO + 0 is less than 1 kcal/mol. Since the vibrational zero 
point energy is ~ 4 kcal/mol, the minimum supports no bound vibrational levels. 
It may be that the flat region of the surface corresponding to the N00 minimum 
is important in photodetachment from N00 — . 

During the course of this work, a large portion of the PES for reaction (2) has 
been mapped out. Since state to state cross sections for reaction (2) are impor- 
tant in the chemistry of high temperature air, these studies will probably be ex- 
tended to permit generation of a new global potential for reaction (2). 

Recently, there have been new experimental studies by Valentini and cowork- 
ers of the product ro-vibrational distributions in the abstraction channel of the 
CH 4 + H reaction. In order to help understand the new experimental results, 
a joint theoretical study is being carried out in collaboration with Muckerman 
( Brookhaven) and Valentini ( Columbia). One feature of the new potential for 
CH 5 , which had not been appreciated before, is that while the saddle point for 
H abstraction is of C 3 t, symmetry ( approaching H collinear with a CH bond), at 
larger H to CH 4 separations the C 3 „ constrained approach is a maximum and the 
energy decreases if the symmetry is lowered. It appears that the minimum energy 
path involves initial approach toward a three fold hole, i.e. collinear with a CH 
bond but on the opposite side of the carbon, with subsequent switch to collinear 
with one CH bond at closer separations. These new features of the potential may 
help explain the correlations of product rotational and vibrational levels oberved 
in the experiments. 

Work is continuing on the HNO system in collaboration with Schatz ( North- 
western). The potential functions for the 1 A', 3 A", and 1 A" surfaces are essen- 
tially complete and preliminary studies have been carried out for the reaction 

OH + N -> H + NO (3) 
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which occurs on the 3 A" surface of HNO. Both the thermal rate constant and 
product ro- vibrational distributions are in good agreement with experiment. These 
studies will be extended to the reaction 


NH + 0 N + OH 

(4a) 

NH + 0 —> H + NO 

(46) 


and will predict the relative probability of reaction (4a) or (4b) (i.e. the product 
branching ratio) and will enable testing of approximate methods for predicting 
product branching ratios. 

In the remaining funding period, it is planned to begin studies of the reactions 
of singlet methylene with N2 and acetylene. The reaction with N2 may be impor- 
tant in “prompt” N0 X and the reaction with acetylene is believed to be impor- 
tant as an initiating step in soot formation. 
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Transport Properties of Gases 
Investigator: Dr. Eugene Levin 


During the time period from 1 January 1993 to 31 August 1993 work continued 
on the determination of the transport properties of gases (e.g., viscosity, diffusion, 
thermal conductivity, etc.) from the collision cross-sections of the constituent 
species. This work was carried out in collaboration with NASA scientists Stall- 
cop, Partridge (and others). 

The results for H-H 2 O completed during the prior reporting period was pre- 
sented at a poster session of the 14th Annual West Coast Theoretical Chem- 
istry/Statistical Mechanics Conference June 17-19, 1993. The paper entitled ”H- 
H 2 0 potential energy surfaces and transport properties” by James R. Stallcop, 
Harry Partridge and Eugene Levin was presented by Dr. Stallcop of NASA. A 
formal submission to the Journal of Chemical Physics is in preparation. Similarly, 
the work on H-H 2 interactions reported during the last time period resulted in a 
new paper, which has been accepted for publication by the Journal of Chemical 
Physics [1]. 

During the present time period, computations were completed for H 2 -H 2 and 
H 2 -N 2 . Cross section calculations for N 2 -N 2 were completed based on a prelimi- 
nary set of potentials. More accurate ab initio potentials have now been obtained 
by Dr. Partridge. Calculation of new cross sections and transport properties are 
awaiting adjustment of the new interaction potentials to form a smooth junction 
with the long-range forms. 

Work has been continuing on the comparison of the semi-classical and quantum 
mechanical methods for the calculation of collision cross sections at low collision 
energies. The results indicate that the difference is insignificant for repulsive po- 
tentials even at low reduced mass of the collision partners and low collision ener- 
gies. However, for bound potentials, the quantum mechanical calculations must 
be used to obtain accurate collision integrals for temperatures below about 1500 
K even for collision partners with a reduced mass as large as 10. This indicates 
that prior work for N-N, 0-0, and N-0 should be revised for such temperature 
regimes. 
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During the next time period the results for the comparison of semi-classical and 
quantum mechanical methods will be documented and submitted for publication. 
The cross section, collision integral and transport property calculations for N2-N2 
will be completed as soon as the adjusted potentials are received from Dr. Stall- 
cop. A long-delayed report on a five-species model for air is still awaiting collision 
cross section or collision integral data for N-e and 0 -e. 

References. 

1 . A b initio potential energy surface for H-H2, H. Partridge, C. W. Bauschlicher, 
Jr., J. R. St allcop and E. Levin, J. Chem. Phys., in press 
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Theoretical characterization of the reaction CH 3 +OH— CH 3 OH — products: 
The 1 CH 2 +H 2 0, H 2 +HCOH, and H 2 +H 2 CO channels 

Stephen P. Walch a) 

ELORET Institute , Palo Alto, California 94303 

(Received 15 September 1992; accepted 10 November 1992) 

The potential energy surface (PES) for the CH 3 OH system has been characterized for the 
! CH 2 +H 2 0, H 2 -l-HCOH, and H 2 + H 2 CO product channels using complete-active-space self- 
consistent-field (CASSCF) gradient calculations to determine the stationary point geometries 
and frequencies followed by CASSCF/intemally contracted configuration-interaction ( CCI ) 
calculations to refine the energetics. The l CH 2 + H 2 0 channel is found to have no barrier. 

The long range interaction is dominated by the dipole— dipole term, which orients the respec- 
tive dipole moments parallel to each other but pointing in opposite directions. At shorter sep- 
arations there is a dative bond structure in which a water lone pair donates into the empty 
a" orbital of CH 2 . Subsequent insertion of CH 2 into an OH bond of water involves a non- 
least-motion pathway. The H 2 + HCOH, and H 2 + H 2 CO pathways have barriers located at 
— 5.2 and 1.7 kcal/mol, respectively, with respect to CH 3 +OH. From comparison of the 
computed energetics of the reactants and products to known thermochemical data it is esti- 
mated that the computed PES is accurate to ±2 kcal/mol. 


I. INTRODUCTION 

The CH 3 +OH reaction has at least six possible prod- 
uct channels: 


CHj+OH + M-CH 3 OH + M, 

(la) 

CH 3 +OH CHjOH + H, 

(lb) 

ch 3 +oh-ch 3 o+h, 

(lc) 

ch 3 +oh-.'ch 2 +h 2 o, 

(Id) 

ch 3 +oh-.hcoh+h 2 , 

(le) 

ch 3 +oh-.h 2 co+h 2 . 

00 


The role of reaction (Id) has been controversial. The room 
temperature rate for the reverse of reaction (Id) has been 
measured by Hatakeyama et aO as ~3xl0 " 12 
cm 3 molecule -1 cm” 1 , while Hack et al. 1 obtained 3.5 
X 10 " 11 cm 3 molecule " 1 cm" 1 , or about an order of mag- 
nitude faster. The latter rate is approximately gas kinetic 
and implies no barrier. Using currently accepted heats of 
formation , 3 and a singlet-triplet splitting for methylene of 
9 kcal/mol, reaction ( Id) is exothermic by 0.7 kcal/mol at 
0 K. Thus the forward reaction is also expected to be very 
fast. 

Dean and Westmoreland 4 have used a variant of Rice- 
Ramsperger-Kassel-Marcus (RRKM) theory called 
QRRK theory to model the product distributions in the 
CH 3 + OH reaction. In this work the parameters for reac- 
tions (lb) and (lc) were based on estimated rates for the 
reverse reaction. The parameters for reaction (If) were 
taken from calculations 5 and the rate for reaction (Id) was 
taken from Ref. 1. According to these studies, at room 
temperature and moderate pressure CH 3 OH is the domi- 


a) Maiting address: NASA Ames Research Center, Moffett Field, 
CA 94035. 


nant product, while at flame temperatures reaction (lb) is 
thought to take over. The HCOH + H 2 channel does not 
appear to have been considered in these studies, though 
theory 5 indicates essentially no barrier with respect to CH 3 
-f-OH. This study indicated that production of *CH 2 [re- 
action 1(d)] is a minor channel. By contrast a model pro- 
posed by Pilling and co-workers 6 which makes use of the 
rate for the reverse of reaction (Id) due to Hack et al. 2 
indicates that reaction (Id) is the dominant channel above 
room temperature. 

Recently Smith 7 has also reported RRKM calculations 
for CH 3 + OH. These calculations as well as the work of 
Pilling et al} have indicated a need for more accurate po- 
tential energy surface (PES) information for the l CH 2 
+ H 2 0, H 2 + HCOH, and H 2 + H 2 CO product channels. 
The most accurate previous theoretical study of these 
channels in the CH 3 OH system was carried out by Hard- 
ing, Schlegal, Krishnan, and Pople 5 using Moller-Plesset 
perturbation theory with a 6-3 1G** basis set. Although 
these calculations were carefully carried out, there is prob- 
ably considerable uncertainty in the energetics; by current 
standards, both the basis set and treatment of electron cor- 
relation can be improved upon. More recently the bond 
dissociation energies of CH 3 OH have been computed by 
Bauschlicher, Langhoff, and Walch 8 using the modified 
coupled-pair functional method. Similar calculations were 
carried out by Pople and co-workers 9 using the G2 
method, which includes some empirical corrections. These 
calculations accurately determined the heats of formation 
of the CH 3 0 and CH 2 OH species, but they did not con- 
sider the portions of the PES leading to the l CH 2 -f-H 2 0 , 
H 2 -bHCOH, and H 2 -bH 2 CO product channels. Thus 
these regions of the PES are reexamined here. 

Qualitative features of the potential energy surfaces are 
discussed in Sec. II, the computational method is discussed 
in Sec. Ill, the results are presented in Sec. IV, and the 
conclusions are given in Sec. V. 


J. Chem. Phys. 98 (4), 15 February 1993 0021 -9606/93 /043163-05$06.00 


© 1993 American Institute of Physics 3163 



Stephen P. Walch: The reaction CH 3 + OH -* CH 3 OH — products 


3165 


TABLE II. Computed energies and zero-point corrections. 


w ch 3 oh^ch 2 o+h 2 

Energy* Zero-point energy b 

CHjOH — 1 15.513 97( —0.54566) 0.052 58 

CH 2 0-H 2 — 1 15.356 03 ( —0.39071) 0.044 07 

CH 2 0+H 2 — 115.472 15( — 0.50189) 0.037 68 

^ CH 3 OH -* *CH 2 4 - H 2 0 


CH 3 OH — 115.513 97 ( —0.54566) 0.052 58 

CH 2 -H 2 0 -115.369 31 (-0.40561) 0.043 66 

CH 2 + H 2 0 -115.359 89 ( —0.39040) 0.038 84 

(c) ch 3 oh^hcoh+h 2 


CHjOH -1 15.505 89( -0.54231) 0.052 58 

HCOH-H 2 -1 15.358 03 ( -0.39640) 0.042 25 

HCOH + H 2 -115.378 23( -0.41465) 0.038 28 




0.0 

91.9 

18.1 


0.0 

82.3 

88.8 


0.0 

85.0 

71.1 


Energy in hartree. The first energy is the CCI energy, while the energy in parenthesis includes a multiref- 
erence Davidson correction and is with respect to - 1 15. hartree 
b Zero-point energy in hartree. 

Relative energy in kcal/mol including zero-point energy and a multireference Davidson’s correction. 


and 'CH 2 + H 2 0 , (denoted as CH 2 0-H 2 , HCOH-H 2 , and 
CH 2 -H 2 0 , respectively) and for the CH 2 + H 2 0 dative 
bonded structure (denoted as CH 2 • H 2 0) are also given in 
Table III. The stationary point corresponding to 
CH 2 • H 2 0 is a minimum on the PES, but there is one very 
small frequency (26 cm -1 ), which corresponds to a hin- 
dered rotation of the CH 2 and H 2 0 with respect to each 
other. 

The computed energy separations discussed here, in 
each case, involve breaking two bonds and forming two 
new bonds; thus the errors in the individual bond strengths 
cancel and the computed energetics are expected to be ac- 
curate. However, in Ref. 2 it was shown that, for calcula- 
tions of about the same quality as reported here, the error 
in the C-O bond strength in CH 3 OH is 6.5 kcal/mol. Thus 
in order to compute energies with respect to CH 3 + OH, 
the experimental 0 K value of 90.2 kcal/mol 18 was used for 
the C-O bond strength. The locations of the H + CH 3 0 
and H + CH 2 OH asymptotes were taken as the best- 
estimate values from Ref. 2. This places H + CH 3 0 and 


TABLE III. Computed saddle point frequencies and rotational constants 

(cm -1 ). 



ch 2 o-h 2 

hcoh-h 2 

ch 2 h 2 o 

ch 2 -h 2 o 

0 )\ 

3195 

4130 

3860 

3723 

"2 

2295 

3199 

3744 

3319 

"3 

1740 

2323 

3187 

3222 

"4 

1574 

1564 

3113 

2091 

"5 

1429 

1498 

1691 

1560 

"6 

1369 

1347 

1515 

1524 


916 

1291 

802 

1119 

"g 

2877/ 

1122 

610 

949 

£t)q 

3278 

934 

337 

703 

"10 

1273 

622 

334 

438 

"11 

1211 

513 

161 

387 

"12 

1065 

1414/ 

26 

1850/ 

A 

3.345 

3.042 

4.117 

4.640 

B 

0.944 

0.839 

0.358 

0.496 

C 

0.863 

0.764 

0.349 

0.478 


H + CH 2 OH at 14.8 and 6.0 kcal/mol above CH 3 + OH, 
respectively. The experimental locations of ‘CHj + HjO 
and CH 2 O-I-H 2 with respect to CH 3 OH were derived from 
the JANAF 3 heats of formation of H 2 , CH 2 0, CH 3 , OH, 
H 2 0(g), and 3 CH 2 , plus a singlet-triplet splitting in CH 2 
of 9.0 kcal/mol and the value for the C-O bond strength in 
CH 3 OH given above. 

From Table I it is seen that the computed 
CH 3 OH— 'CHj + HjO separation is 0.7 kcal/mol 
smaller than experiment, while the computed 
CH 3 0H — CH 2 0 +H 2 separation is 0.6 kcal/mol smaller 
than experiment. It is also seen that the computed results 
of Harding et al . 5 are 5.4 kcal/mol larger and 2.3 kcal/mol 
smaller, respectively, for the same separations. The com- 



FIG. 1. Schematic diagram of the potential energy surface for CH 3 + OH. 
The location of the CH 3 + OH asymptote with respect to CH 3 OH is taken 
from experiment, while the locations of the CH 2 OH + H aand CH 3 0 + H 
asymptotes are from previous calculations. 
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TABLE V. Computed dipole and quadrapole moments. 3 




Q*x 

Q» 

Qn 

h 2 o 

0.768 

-1.987 

1.955 

0.032 

ch 2 '^. 

0.679 

1.187 

0.496 

-1.683 


"Properties are in a.u. Quadrapole moment is with respect to the center of 
mass. The molecule is in the YZ plane with the C 2 axis in the Z direction. 


equilibrium geometries and oriented with the planes of the 
molecule parallel to each other and the CO bond perpen- 
dicular to both molecular planes in the orientation shown 
in Fig. 2. For this geometric orientation, ^ was varied 
and the resulting energies at the CCI level are also included 
in Fig. 2. 

The main features of Fig. 2 are a shallow minimum 
followed by a small barrier to formation of CH 3 OH. How- 
ever, the barrier is below the ! CH 2 -(~H 20 asymptote and 
therefore the bottleneck on the vibrationally adiabatic 
curve is expected to occur in the entrance channel region. 
The main feature responsible for the entrance channel bot- 
tleneck is the building in of bending modes, which arise 
from electrostatic (dipole-dipole, dipole-quadrapole, and 
quadrapole-quadrapole) interactions. In order to define 
this interaction, dipole and quadrapole moments (about 
the center of mass) were computed for J CH 2 and H 2 0 and 
are given in Table V. The experimental values 19 for H 2 0 
(in a.u.) are jx = 0.73, Q xx = — 1 859, Q yy — — 1955, and 
Q zz = —0.097; which are in reasonable agreement with the 
computed values. 

V. CONCLUSIONS 

The potential energy surface (PES) for the CH 3 OH 
system has been characterized for the 1 CH 2 -hH 2 0, H 2 
+ HCOH, and H 2 +H 2 CO product channels using 
complete-active-space self-consistent field (CASSCF) gra- 
dient calculations to determine the stationary point geom- 
etries and frequencies followed by CASSCF/intemally 
contracted configuration-interaction (CCI) calculations to 
refine the energetics. 

The H 2 -f-H 2 CO, and H 2 -{-HCOH pathways have bar- 
riers located at 1.7 and —5.2 kcal/mol with respect to 
CH 3 + OH. The l CH 2 +H 2 0 channel is found to have no 
barrier in the absence of vibrational zero-point effects. 
However, the long range interaction is dominated by a 


dipole-dipole term and the zero-point effects due to this 
interaction are expected to lead to a bottleneck on the 
vibrationally adiabatic minimum energy path. The ! CH 2 
+ H 2 0 asymptote is computed to be 1.4 kcal/mol below 
CH 3 + OH. Thus all three of these channels are expected to 
be accessible at moderate temperatures. 

From comparison of the computed energetics of the 
reactants and products to known thermochemistry it is 
estimated that the computed PES is accurate to ± 2 kcal/ 
mol. 
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Computed barrier heights for H+CH 2 O^CH 3 O^CH 2 OH 

Stephen P. Walch a) 

ELORET Institute, Palo Alto t California 94303 

(Received 4 September 1992; accepted 28 October 1992) 

The barrier heights (including zero-point effects) for H + CH 2 0 — CH 3 0 and 
CH 3 0-»CH 2 0H have been computed using complete active space self-consistent field 
(CASSCF) /gradient calculations to define the stationary point geometries and harmonic 
frequencies and internally contracted configuration-interaction (CCI) to refine the energetics. 
The computed barrier heights are 5.6 and 30.1 kcal/mol, respectively. The former barrier 
height compares favorably to an experimental activation energy of 5.2 kcal/mol. 


The CH 3 0 radical is important in combustion and at- 
mospheric chemistry. 1 In addition, recent experiments 
have studied individual rovibrational states of CH 3 0 above 
the H-CH 2 0 dissociation threshold by stimulated emission 
pumping. 2 The interpretation of these experiments depends 
critically on the barriers to the decomposition of CH 3 0 to 
H + CH 2 0 and rearrangment to CH 2 OH. Previous studies 
of these barrier heights include the work of Saettff, 
Radom, and Schaefer, 3 who used Mdller-Plesset perturba- 
tion theory through third order with small basis sets up 
through 6-3 1G**. More recently, Page, Lin, He, and 
Choudhury 4 reported somewhat more accurate calcula- 
tions using multireference configuration-interaction 
(MRCI) with polarized triple zeta basis sets. In the 
present paper more accurate calculations are reported 
which make use of larger basis sets and extensive Cl. 

Two different basis sets were used in this work. For the 
complete active space self-consistent field (CASSCF) gra- 
dient calculations the polarized double zeta set of Dunning 
and Hay 5 was used. The basis set for C and O is a (9s5p)/ 
[3s2/>] basis augmented by a single set of 3 d functions with 
exponents of 0.75 and 0.85 for C and O, respectively. The 
H basis is (4s) /[2s] augmented with a single set of 2 p 
functions with exponent 1.00. The basis set used in the Cl 
calculations is the Dunning correlation consistent triple 
zeta double polarization basis set. 6 This basis is [4s3pld\J] 
for C and O and [3s2p \d] for H and is described in detail 
in Ref. 6. 

The calculations were carried out in C s symmetry for a 
wave function of 2 A' symmetry. From the Hessian matrix 
in the diagonal representation, the stationary points ob- 
tained in C s symmetry are the minimum for CH 2 0-f H, the 
saddlepoint for CH 2 0 + H — CH 3 0, the CH 3 0 minimum, 
and the saddlepoint for CH 3 0->CH 2 0H. For CH 3 0 the 
2 A f state corresponds to one component of the 2 E state in 
Qu symmetry. According to the calculations of Saebd et 
aly 3 the Jahn-Teller stabilization energy in this system is 
0.56 kcal/mol and the Jahn-Teller splitting is 0.12 kcal/ 
mol, with the 2 A’ component lower. CH 2 OH is found to 
have no symmetry and the use of a mirror plane of sym- 
metry would not be appropriate. However, the separation 
between CH 3 0 and CH 2 OH has been accurately computed 


a) Mailing address: NASA Ames Research Center, Moffett Field, Califor- 
nia 94035. 


by Bauschlicher, Langhoff, and Walch 7 and further calcu- 
lations for CH 2 OH were not carried out here. 

The CASSCF calculations consisted of five electrons 
and five orbitals. The active electrons for CH 3 0 included 
the CO bond pair, one CH bond pair (the CH bond in the 
molecular plane), and the a' O 2 p like orbital. These or- 
bitals correspond to the CO a and v bonds plus a H Is 
orbital for CH 2 0 + H. The remaining electrons, which are 
inactive, include the other two CH bond pairs, the a” O 2 p 
like orbital and the O Is, O 2s, and C Is like orbitals. In 
generating the set of reference configurations for the sub- 
sequent internally contracted Cl (CCI) calculations, no 
more than two electrons were permitted in the weakly oc- 
cupied CASSCF orbitals. All but the O Is and C Is elec- 
trons were correlated in the CCI calculation. A multiref- 
erence analog of Davidson’s correction 8 was added to the 
CCI energies. 

The CASSCF/gradient calculations used the sirius/ 
ABACUS system of programs, 9 while the CCI calculations 
were carried out with MOLPRO. 10,11 The calculations were 
carried out on the NASA Ames Cray Y-MP. 


TABLE I. Computed stationary point harmonic frequencies (cm '), 
rotational constants (cm 1 ), bond lengths (A), and bond angles (de- 
grees). 



CHjO 

CH,O^CH 2 OH 

h-ch 2 o 

h + ch 2 o 


3216 

3261 

3190 

3174 


2921 

2378 

1657 

1766 


1608 

1556 

1471 

1589 

<y 4 

1471 

1101 

1220 

1213 

“>5 

1103 

1037 

605 



997 

2254/ 

1435/ 


Ci)j 

3284 

3391 

3297 

3272 


1491 

1159 

1281 

1304 

£l)q 

1050 

1015 

770 


A 

5.180 

5.718 

3.927 

9.548 

B 

0.885 

0.929 

1.007 

1.263 

C 

0.884 

0.900 

0.962 

1.115 


1.092 

1.086 

1.094 

1.097 

r C-H h 

1.122 

1.275 

1.668 


'c-o 

1.412 

1.418 

1.259 

1.226 

'b-H* 


1.237 

2.292 


Z HyCHy 

111.5 

119.1 

117.1 

117.2 

/ H.CH* 

108.2 

117.1 

91.2 


Z H„CO 

111.7 

116.5 

120.3 

121.4 

ZH,CO 

105.3 

54.4 

102.2 
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Cl (ICCI) methods. The double, -quartet crossing region ts found to h- T/ SCF<CASSCF > and '"'anally contracted 
and from products by a barrier of % 23 keal/mof 1 separated from reactants by a barrier of *18 keal/mol 


1. Introduction 


The reaction 


CH(X 2 n) + N 2 (X 'S g + ) 
-HCN(X l I + )-|-N( 4 S) 


0 ) 


? bclicv( ; d t0 be the initiating step for formation of 

prompt” NO in hydrocarbon combustion [ 11 . 
Measurements of the rate of reaction ( 1 ) have been 
summarized by Miller and Bowman [ 1 ]. There is a 
considerable variation of the measured rate among 
the different experiments; from the rate expressions 
the activation energy varies from 1 1 keal/mol to 
*20 keal/mol. The A factor for the reaction is re- 
ported f 1 ] to be about a factor of 100 smaller than 
would be expected for a similar reaction which was 
spin allowed. 

Manaa and Yarkony [2,3] have used a con- 
strained analytical gradient search algorithm to lo- 
cate regions of the doublet-quartet crossing hyper- 
surface and have also computed the spin-orbit matrix 
element coupling the two surfaces. Thev find that the 
lowest energy point on the doublet-quartet crossing 
>persurface corresponds approximately to a Co, 
configuration (CH bond perpendicular to the mid- 

; A “ 6 a L ^ CSS; NASA Research Center, Moffett Field, 
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point of a stretched N 2 bond). This region is esti- 
mated to be ^ 7.5 keal/mol above the CH ( 2 n ) + N, 
reactants. These workers also characterized a C 2v " 
minimum on the doublet surface, denoted min(C-> ) 
which is similar to the structure of the crossing rc- 
gron except for a shorter NN bond, and another min- 
imum, denoted as min (dative), on the same surface 
These authors suggest the formation of an interme- 
diate metastable complex structurally similar to 
min(C 2v ) which repeatedly traverses the doublet- 
quartet crossing region resulting ultimately in the 
formation of HCN + N( 4 S) products. Fig. 1 of ref. 
[3] suggests a pathway in which the barrier to re- 
action is at the doublet-quartet crossing, but no bar- 
riers separate this region from reactants or products. 
Bair [4], in an unpublished account, suggested the 
formation of min (dative) involved no barrier but 
no detailed calculations had been carried out to de- 
termine whether entrance or exit channel barriers 
existed. 

More recently, Martin and Taylor [5] have ap- 
plied the same methods used in the present work to 
this system. They considered an entrance channel 
pathway in which CH( 2 I7) + N 2 initially form 
min(dative) and then convert to min(C 2v ) leading 
to the doublet-quartet crossing region studied bv 
Manaa and Yarkony. However, they find a verv large 
barrier (> 40 keal/mol) separating min (dative) 
from mm(C 2v ). They also find an exit channel bar- 

00 © 1993 Elsevier Science Publishers B.V. All rights reserved. 
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r.er of the same magnitude. They conclude “that the 
reaction mechanism is different from what has been 
previously suggested” (i.e. Manaa and Yarkony sig- 
nificantly underestimated the barrier between 
min (dative) and min(C 2v )) ; however, they were not 
able to find an alternative low-energy pathway 
The failure to find a low-energy pathway for in- 
”“" of CH <’ n > i»'° Nr is surprising, since 
LH( 1 1 ' ls Slm, lar electronically to 'CH 2 , which typ- 
ically inserts into bonds with little or no barrier. 
However, previous experience with the CH( 2 rj) + H 2 
reaction [6,7] shows that insertion by a least-motion 
path involves a high barrier, while non-least-motion 
pathways involve a small barrier. The pathway in- 
volved in the latter reaction may be understood in 

of the orbltal P hase continuity principle 
(OPCP) proposed by Goddard [8], In section 2 we 
show that use of the OPCP predicts that the CH ad- 
dition pathway studied by Martin and Taylor (trans 
orientation) will have a high barrier, while a cis ori- 
entation leads to addition with a low barrier. 


2. Qualitative features 

The X 2 n state of CH may be represented by 

( 2 ) 



i.e. the two orbitals corresponding to the C lone pair 
are singlet paired. If these orbitals are solved for self 
consistently in a generalized valence bond wave- 
function the overlap integral between them is *0.7. 
Thus, the X 2 U state of CH may be characterized as 
a singlet biradical, but the orbitals of the lone pair 
have a substantial overlap, which must be main- 
tained while inserting into a single bond if the pro- 
cess is to occur without a substantial barrier. In 
terms of multiconfigurational self-consistent-field 
(MCSCF) theory, the biradical character in the X 2 n 
state of CH arises because of a near degeneracy effect 
between an sp hybrid lone pair and an empty 2p-like 
orbital. 

A complete discussion of the OPCP is beyond the 
scope of this Letter and the reader is referred to the 
original paper by Goddard [8] and also to the dis- 


cussion of the CH 2 n + H 2 reaction in ref. [7], In the 
case of the CH 2 n + N 2 reaction, for approach in a 
trans orientation (the pathway considered by Mar- 
tin and Taylor) (eq. (3), see below), it is seen that 
transfer of the C lone pair to become a CN bond pair 
leads to a nodal plane between the upper N 2p or- 
bital of the N 2 in plane tt bond and the other N 2p 
orbital, which evolves into a singly occupied orbital 
on C. (In eqs. (3) and (4) below positive orbital 
obes are indicated by solid lines, negative orbital 
lobes are indicated by dashed lines, and the straight 
dashed line in eq. (3) indicates a nodal surface which 
results from overlap of orbitals of opposite phase.) 
Thus, this pathway involves breaking of one bond 
and is expected to have a high barrier, as Martin and 
Taylor found. The alternative cis approach path (eq. 
(4)), permits orbital phases consistent with trans- 
fenng both bond pairs simultaneously, and this path- 
way should involve a much smaller barrier, 



3. Computational details 


Two different basis sets were used in this work. 
For the CASSCF geometry optimizations, which were 
performed using analytical first and second deriva- 
tive methods, the polarized double zeta set of Dun- 
ning and Hay [9] was used. The basis set for C and 
N is a (9s5p)/[3s2p] basis augmented by a single 
set of 3d functions with exponents of 0.75 and 0.80 
for C and N, respectively. The H basis is (4s)/ [2s] 
augmented with a single set of 2p functions with ex- 
ponent 1 .00. The basis set used in the Cl calculations 
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Table 2 


Computed stationary point energies a> 


ICCI(ICCI + Q+ 147.) ZPE b| 1CCI + Q+ 147. +ZPE 

reac — 147.58272( —0.59243 ) 0.01119 

dspl -147.54932( -0.56830) 0.01603 

dm in 1 — 147.57896( —0.59738) 0.02014 

dsp2 — 147.56897( —0.59201 ) 0.01768 

dmin2 - 147.60531 ( -0.62126) 0.01976 

qminl - 147.57758( -0.59285) 0.01843 

<)SPI — 147. 54642(— 0.56127) 0.01643 

a> Energies in E h unless otherwise noted. 
b> Zero-point energy based on harmonic CASSCF frequencies. 
c) Relative energy (kcal/mol). 


-0.58124 

0.0 

-0.55227 

18.2 

-0.57724 

2.5 

-0.57433 

4.3 

-0.60150 

-12.7 

-0.57442 

4.3 

-0.54484 

22.8 




Fig. 2. Schematic representation of the minimum energy path for 
reaction ( I ). The energies are from 1CCI calculations correlating 
eleven electrons. The location of the doublet— quartet surface 
crossing is taken from the work of Manaa and Yarkony [2,3). 


Fig. 1. Stationary point geometries along the minimum energy 
path for reaction ( 1 ). Note that dspl, dminl, and dsp2 are non- 
planar (see table I ). 


imposed in the calculations. 

From fig. 1 it is seen that the CH and N 2 approach 
each other in a cis orientation consistent with eq. (4). 
An unexpected feature of this surface is the presence 
of the minimum denoted by dminl. This is a short 
NN structure with C s symmetry' (the CH bond lies 
in a plane perpendicular to the C-N,-N 2 plane and 
passing through the C and midpoint of the NN 
bond). From fig. 2 it is seen that dminl is only a 


shallow minimum connected by dsp2 to dmin2, 
which is a planar C 2v structure with a long NN bond. 
dmin2 is the same as Manaa and Yarkony’s min ( C 2v ) 
and Martin and Taylor’s D3 structure. Thus, the key 
feature of the present work is a pathway connecting 
CH 2 n + N 2 to Manaa and Yarkony’s min(C 2v ) with 
a barrier of only %I8 kcal/mol. By contrast, the 
pathway via Manaa and Yarkony’s min (dative) 
structure involves a trans approach and has a barrier 
> 40 kcal/mol in the calculations of Martin and Tay- 
lor. It is probable that, while formation of the 
min (dative) structure occurs with no barrier, the 
large barrier separating it from min(C 2v ) means that 
this pathway does not lead to products. 


217 



Theoretical Characterization of the Reaction 


NH 2 + NO — * products 


Stephen P. Walch a 
ELORET Institute 
3788 Fabian Way 
Palo Alto, Ca. 94303 

Abstract. The potential energy surface for NH 2 + NO has been characterized 
using complete active space self consistent field (CASSCF)/ derivative calculations 
to determine the stationary point geometries and frequencies followed by internally 
contracted configuration interaction (ICCI) calculations to determine the energetics. 
Production of N 2 + H 2 0 is found to involve a complex mechanism. The initially 
formed NH 2 NO undergoes a 1,3-hydrogen shift to give an HNNOH isomer ( with 
the substituents trans about the NN bond and cis about the NO bond) which 
undergoes subsequent cis-trans isomerizations about the NN and NO bonds before 
decomposing to N 2 -f- H 2 0. The saddle point for production of N 2 + II 2 0 has an 
approximately rectangular arrangement of one H atom, the two N atoms, and the O 
atom. This process does not involve a barrier with respect to NH 2 + NO. Formation 
of HN 2 + OH can occur from any of the isomers of HNNOH with no barrier, but 
the overall process is endothermic by 0.7 kcal/mol (based on the computed A H^ 
(OK) of HN 2 ). 

The results obtained in this work are qualitatively the same as previous work, 
but both the stationary point geometries and energies should be more reliable due 
to the use of larger basis sets and more extensive inclusion of electron correlation 
effects. 

“Mailing Address: NASA Ames Research Center, Moffett Field, CA 94035. 
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I. Introduction. 


The thermal deNO x process [1-4], in which the addition of NH3 to combustion 
processes reduces the production of NO x , may be useful in reducing N0 Z emissions 
from jet engines. The reaction of NH2 + NO is believed to be a key reaction in this 
process, but the product branching ratio at flame temperatures is not known. The 
two major channels which have been considered are: 


nh 2 + no —> n 2 + h 2 o 

(1) 

NH 2 + NO — > HN 2 + OH 

(2) 


The experimental literature on this reaction has been discussed by Silver and Kolb 
[5]. The product branching ratio a = k2/(ki + k2) (where ki and k2 are rate co- 
efficients for reactions (1) and (2), respectively) has ranged from > 0.65 to < 0.13 
at room temperature [6-10], although most recent determinations [9-10] have con- 
verged toward lower values (0.13 [9], < 0.13 [10], and 0.12 [5]). Note that in ref. 
5, Silver and Kolb reevaluated their earlier results [6] to correct for OH production 
from the reaction 

NH + NO -> N 2 + OH (3) 

The value of the branching ratio (a) is critical in the mechanism of the thermal 
deNO x process, because production of radical species via reaction (2) feeds chain 
reactions. Attempts to model the thermal deNO x process have required a value of a 
of 0.29-0.4 to model measurements under combustion conditions [11-16]. The work 
of Kimball- Linne and Hanson [17] implied that a varies with temperature from 0.48 
at 1050K to >0.8 at 1400K. No direct experimental determinations exist at flame 
temperatures. 
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It is also possible to obtain the brandling ratio ( a ) from theory. This would 
require accurate characterization of the potential energy surface (PES) for NH2 + 
NO. There have been several previous theoretical studies of the PES for NH2 + 
NO. Casewit and Goddard [18] used the generalized valence bond - configuration 
interaction (GVB-CI) method with a polarized double-zeta basis set to study nine 
isomers of formula N2H2O. In these studies the geometries were obtained using 
restricted Hartree-Fock (RHF) theory with a 4-31G basis set. One minor problem 
in this work was that a planar structure was found for NH2NO (N-nitrosamide). 
Later studies by Harrison et al. [19] showed that addition of polarization functions 
leads to a non-planar structure. 

Melius and Binkley [20] carried out an extensive series of calculations on this 
system using the bond additivity corrected / fourth order M0ller-Plesset perturba- 
tion theory ( BAC-MP4 ) method with a split- valence plus polarization basis set. 
In addition to the large bond additivity corrections (10-20 kcal/mol per bond) ap- 
plied in these calculations, there are also corrections for spin contamination in the 
unrestricted Hartree-Fock (UHF) calculations. In previous studies on the PES for 
NH + NO [21] it has been shown that, while the BAC-MP4 method gives results 
in qualitative accord with the results of multireference configuration interaction 
calculations, there are differences in the detailed energetics as large as 8 kcal/mol. 

In a subsequent paper, Harrison and coworkers [22] carried out a study using 
MP2 through MP4 theory with a 6-31G* basis set. In addition to characterizing 
the minima on the surface, they also located a number of saddle points. These cal- 
culations differed from those of Melius and Binkley in that no empirical corrections 
were used. Comparison of the results obtained using this approach to MRCI [23] 
results for the singlet surface of CH3OH indicates errors as large as 5 kcal/mol in 
the detailed energetics. 
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In order to improve on the energetics of the NH2 + NO reaction calculations 
were carried out here using CASSCF/derivative and ICCI methods with large basis 
sets. Section II presents some qualitative features of the calculations, Section III 
discusses the technical details of the calculations, Section IV presents the results 
and the conclusions are given in Section V. 

II. Qualitative Features. 

The ground state of NH2 is 2 Bi, while the 2 Ai state is a low-lying excited state. 



Combining the ground state of NH2 with the 2 II state of NO is most favorable 
for a non-planar approach 


leading to the NH2NO species (denoted as mini) with a pyramidal geometry 
about the NH 2 N atom. 

Fig. 1 shows the electronic structure of the stable minima on the potential energy 
surface. In order to consistently describe the making and breaking of bonds over the 
whole surface, it is necessary to correlate the two NH bonds, the NN a bond, the 
NO (7 bond, and the NO 7 r bond of mini. Note that these orbitals ultimately evolve 
to an NN a bond, two NN 7 r bonds, and two OH a bonds in the N2 + H2O product 
channel. This choice of correlated orbitals is also able to describe dissociation to 

HN 2 + OH. 
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While mini is non-planar, the remaining structures shown in Figl as well as the 
saddle points connecting them were found to be planar. 

III. Computational Details. 

Two different basis sets were used in this work. For the CASSCF derivative 
calculations the polarized double zeta set of Dunning and Hay [24] was used. The 
basis set for N and 0 is a (9s5p)/[3s2p] basis augmented by a single set of 3d 
functions with exponents of 0.80 and 0.85 for N and 0, respectively. The H basis 
is (4s)/ [2s] augmented with a single set of 2p functions with exponent 1.00. The 
basis set used in the Cl calculations is the Dunning correlation consistent triple zeta 
double polarization atomic natural orbital basis set [25]. This basis is [4s3p2dlf] 
for N and 0 and [3s2pld] for H and is described in detail in Ref. 25. 

For the planar structures ( min2-min4) the CASSCF /derivative calculations had 
10 electrons distributed among 10 orbitals (denoted by 10/10). The active electrons 
correspond to the electrons in singly occupied orbitals in Fig. 1. Calculations were 
also carried out using a smaller CASSCF calculation which had 8 electrons dis- 
tributed among 6 orbitals (denoted by 8/6). For mini the active electrons include 
the upper NH <x bond, the in plane O 2p lone pair, the NO 7 r bond, and the N 2p 
like lone pair. This active space is found to be adequate for describing structures 
mini through min4 and the saddle points connecting them, but cannot describe dis- 
sociation to the N 2 + H 2 O or HN 2 + OH product channels. Finally, a calculation 
was carried out to determine the NN bond strength for the mini structure. The 
CASSCF calculation in this case had 6 electrons distributed among 6 orbitals (de- 
noted by 6/6). The active electrons correspond to the NN a bond, the NO a bond, 
and the NO 7 r bond for mini. In the subsequent internally contracted configuration 
interaction (ICCI) calculations, the N Is and O Is electrons were not correlated. 
These Cl calculations were based on the 10/10, 8/6, or 6/6 CASSCF active spaces, 
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respectively. The reference space for the IGCI calculations was restricted to no more 
than two electrons in the weakly occupied CASSCF orbitals. 

The CASSCF /derivative calculations used the SIRIUS/ABACUS system of pro- 
grams [26], while the ICCI calculations were carried out with MOLPRO [27-30]. 
A multi-reference analog of the Davidson’s correction [31] was added to the ICCI 
energies and is denoted by +Q. 

IV, Discussion 

Tables I, II, and III give computed ICCI energies, zero-point energies (based 
on harmonic CASSCF frequencies) and relative energies in kcal/mol based on the 
(8/6), (10/10), and (6/6) CAS active spaces, respectively. The energetics given in 
Tables I-III are also shown graphically in Fig. 2, while the geometries corresponding 
to the stationary points are given in Fig. 3. 

The relative energies in Table I are given both with respect to min2 and with 
respect to NH 2 + NO. The location of the NH 2 + NO asymptote is derived from 
Table II using the computed separation between min2 and H 2 + N 2 0 plus the 
experimental separation [32] between H 2 + N 2 0 and NH 2 + NO. Since the H 2 + 
N 2 0 and min2 structures involve the same number of bonds ( sometimes referred 
to as an isodesmic reaction), this computed separation is more accurate than that 
obtained by direct computation of the energy to dissociate mini to NH 2 + NO, as 
in Table III. Table II also gives the energy of the H 2 0 + N 2 asymptote, computed 
with respect to H 2 + N 2 0. The computed separation of 124.0 kcal/mol compares 
well to an experimental separation [32] of 124.5 kcal/mol. Comparing the energy for 
mini given in Table I to the directly computed value given in Table III, indicates 
the directly computed NN bond strength in the mini structure is too small by 
5.9 kcal/mol. This is a reasonable value, since an estimate of 7.2 kcal/mol is obtained 
from 2/3 of the error in the N 2 dissociation energy with the same basis set. 
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From Tables I and II the separation between min2 and min3 is 1.5 kcal/mol with 
the 8/6 CASSCF active space and 1.2 kcal/mol with the 10/10 CASSCF active 
space. This indicates that the 8/6 active space provides an adequate description of 
the structures mini to min4, as asserted in Section III. 

From Fig. 2 it is seen that the saddle points in the pathway leading to N 2 + 
H 2 O are all lower in energy than the NH 2 + NO asymptotic energy. The largest 
barrier is for the saddle point separating min3 and min4. From Fig. 3 it is seen that 
this is a planar saddle point, and may be thought of as an inversion process about 
Nl. An alternative pathway involving rotation about the NN bond is found to be 
much higher in energy, in agreement with Fig. 1, where it is seen that min2 and 
min3 both have an NN double bond. The saddle point separating min2 and min3 
corresponds to rotation about the NO bond and has a low barrier, in agreement 
with the bonding structures given in Fig. 1. Saddle point spl for the 1,3-hydrogen 
transfer also has a large barrier, but smaller than for sp3. From Fig. 1 it is seen 
that this process is complex, involving the simultaneous transfer of an NH bond 
pair to become an OH bond pair and transfer of an 0 2p like lone pair to become 
a lone pair on Nl. 

From Fig. 2 it is also seen that structures min2, min3, or min4 can dissociate 
to HN 2 + OH with no barrier other than the exothermicity. From the computed 
A HJ(OK) of HN 2 [33] and the known heats of formation of the other species [32] 
reaction (2) is computed to be endothermic by ~ 1 kcal/mol. 

Reaction (1), on the other hand is exothermic by 124.5 kcal/mol and has no 
barriers with respect to NH 2 + NO. This would suggest that reaction (1) should 
dominate at low temperatures, but reaction (2) may also be important at higher 
temperatures. This is consistent with the bulk of the recent experimental infor- 
mation, particularly ref. 17 which implies that a is 0.48 at 1050K and > 0.8 at 
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1400K. 


The geometric parameters of the structures shown in Fig. 3 are given in Table IV. 
Comparing the geometries for the minima to the corresponding structures found by 
Harrison et al. [22] shows the expected trends, i.e. the CASSCF bond lengths are 
consistently larger than the SCF bond lengths. The saddle point structures found 
here differ more from the structures found by Harrison et al. This probably reflects 
the expected difficulty in describing the process of elongating chemical bonds with 
SCF theory. 

Table V gives computed frequencies and rotational constants for the stationary 
points obtained in this work. Comparing the frequencies obtained here to those ob- 
tained by Harrison et al. indicates the expected trends, i.e. the CASSCF frequencies 
are all smaller than the SCF frequencies. Once again this well known effect derives 
from incorrect dissociation of chemical bonds at the SCF level of theory. This is 
commonly corrected for by an empirical scaling of the frequencies by 0.89. This is 
found to be unnecessary for CASSCF frequencies [21], 

Table VI compares the energetics obtain in the present work with those obtained 
by Melius et al. [20] and by Harrison et al. [22]. In Table VI the results obtained 
here are referenced to NH 2 + NO. This is done in two ways. In the results denoted 
by refl the energies were computed with respect to H 2 + N 2 0 and the NH 2 + 
NO asymptote was positioned from experiment. In the results denoted by ref2 
the energies are directly computed with respect to NH 2 + NO. As discussed above 
the results denoted refl are more reliable. The ref2 results are also included for 
comparison to the work of Harrison et al., which was done in the same way. 

Comparing our refl results to the BAC-MP4 results of Melius et al., it is seen that 
the BAC-MP4 results for the stationary points mini - sp4 vary from 5.3 kcal/mol 
deeper (spl) to -0.3 kcal/mol shallower (sp3). Since the BAC-MP4 method uses 
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rather large (10-20 kcal/mol) bond additivity corrections, it is encouraging that 
the present results are on the average in fairly good agreement with the BAC-MP4 
results ( average error 2.4 kcal/mol more shallow). It is a bit surprising that the 
relative energies of structures mini - sp4 are not in better agreement with the 
present results. What would be expected is that the relative energies of structures 
mini - sp4 should be quite accurate in the present calculations but might be shifted 
somewhat with respect to NH 2 + NO. It is probable that the energy with respect 
to NH 2 + NO is accurate to ± 2 kcal/mol and the relative energies should be 
good to similar accuracy. This suggests that there are disagreements between the 
present work and the BAC-MP4 results that are larger than the expected error in 
the present calculations. On the other hand the BAC-MP4 results do give a realistic 
semiquantitative picture of the PES for this system. 

We now compare to the work of Harrison et al. The column of Table VI la- 
beled ICCI+Q is computed from direct calculation with respect to NH 2 + NO and 
without zero-point energy for comparison to the results quoted in their paper. Com- 
paring these numbers to those quoted by Harrison et al. it is seen that our results 
are deeper by 3.4 - 10.6 kcal/mol ( average 7.0 kcal/mol). Note that our results 
obtained by direct calculation with respect to NH 2 + NO are in turn shallower by 
5.9 kcal/mol than the results based on a calculation with respect to H 2 + N 2 0 and 
use of the experimental separation between H 2 + N 2 0 and NH 2 + NO. Thus, our 
results are deeper than those of Harrison et al. by ^ 13 kcal/mol on an average 
basis. This difference reflects both the use of a larger basis set and a more con- 
sistent CASSCF/ICCI treatment of the electron correlation problem in the present 
calculations. 

Finally, we compare to the work of Casewit and Goddard [18]. Casewit and 
Goddard used a GVB-CI with the Dunning polarized valence double zeta basis set 
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[24], Their heats of formation were computed with respect to the sum of separate 
calculations on N 2 and H 2 O. Since this is an isodesmic approach, their results should 
be compared to our results computed with respect to H 2 + N 2 O. From Table VI it is 
seen that there are substantial errors, ranging from 13.8 kcal/mol to 15.1 kcal/mol 
for mini - min4. These errors are surprisingly large. Part of the error must be due 
to the relatively small basis set, but another possible problem is the use of separate 
calculations on N 2 and H 2 0. A supermolecule calculation would have been more 
appropriate and might have given more accurate heats of formation. 

V. Conclusions. 

The PES for NH 2 + NO has been characterized using CASSCF/derivative meth- 
ods to locate the stationary points followed by CASSCF/ICCI calculations using 
a [4s3p2dlf/3s2pld] correlation consistent atomic natural orbital basis set to de- 
termine the energetics. Production of N 2 + H 2 O is found to involve a complex 
mechanism, which, however, has no barrier with respect to NH 2 + NO. This path- 
way is exothermic by 124.5 kcal/mol. Production of HN 2 + OH can occur with no 
barrier other than the exothermicity from any of three isomers of HNNOH. Using 
the computed heat of formation of HN 2 this channel is endothermic by ~ 1 kcal/mol. 
This suggests that the N 2 + H 2 O channel will dominate at low temperatures, but 
the HN 2 + OH channel will be important at higher temperatures, in agreement 
with the bulk of recent experimental information. 

The present calculations are in semiquantitative agreement with the BAC-MP4 
results of Melius and Binkley, but there are differences which are larger than the 
estimated uncertainty of the present results. The present PES is deeper with respect 
to NH 2 + NO by ~ 13 kcal/mol as compared to the results of Harrison et al. This 
difference arises from two causes. First the present calculations use larger basis 
sets and more extensive treatment of electron correlation. Secondly it is found that 
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computing the location of the NH2 + NO asymptote by computing the energy of 
the bound species with respect to H2 + N2O and then using the experimental H 2 
+ N2O to NH2 + NO separation is more accurate than direct computation with 
respect to NH2 + NO as was done in the Collins et al. work. The latter method is 
found to lead to an error of & 6 kcal/mol ( NH2 + NO too low with respect to the 
rest of the surface). 

Acknowledgement. SPW was supported by NASA cooperative agreement num- 
ber NCC 2 - 478 . 
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Table I. Computed energies and zero-point corrections. 

ICCI calculations based on an (8/6) CASSCF 



Energy® 

zero-point energy*’ 

A E c 



mini 

-185.56699 

0.03494 

-0.53205 

1.7 

-44.0 

spl 

-185.51535 

0.03040 

-0.48495 

31.3 

-14.4 

min2 

-185.57044 

0.03562 

-0.53482 

0.0 

-45.7 

sp2 

-185.55371 



« 10.0 

« -35.7 

min 3 

-185.56798 

0.03554 

-0.53244 

1.5 

-44.2 


a ICCI +Q energy in Eh- 

k zero-point energy in Eh> obtained from the harmonic CASSCF frequencies. 
c First column ICCI +Q energy plus zero-point energy + 185. Second column 
energy in kcal/mol relative to min2. Third column energy in kcal/mol relative to 
NH 2 + NO. 
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Table II. Computed energies and zero-point corrections. 

ICCI calculations based on an (10/10) CASSCF 



Energy® 

zero-point energy 6 

A E c 



h 2 + N 2 0 

-185.57072 

0.02075 

-0.54997 

0.0 

-47.2 

h 2 o + n 2 

-185.69880 

0.02643 

-0.67237 

-76.8 

-124.0 

min2 

-185.58060 

0.03303 

-0.54757 

1.5 

-45.7 

min3 

-185.57839 

0.03278 

-0.54561 

2.7 

-44.5 

sp3 

-185.51521 

0.02873 

-0.48648 

39.8 

-7.4 

min4 

-185.57971 

0.03242 

-0.54729 

1.7 

-45.5 

sp4 

-185.53541 

0.02568 

-0.50973 

25.3 

-21.9 


“ ICCI + Q energy in Eh- 

^ zero-point energy in Eh> obtained from the harmonic CASSCF frequencies. 
c First column ICCI +Q energy plus zero-point energy + 185. Second column 
energy in kcal/mol relative to H 2 + N 2 O. Third column energy in kcal/mol relative 

to NH 2 + NO. 
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Table III. Computed energies and zero-point corrections. 


ICCI calculations based on an (6/6) CASSCF 



Energy* 

zero-point energy 6 

A E c 


NH 2 + NO 

-185.50010 

0.02278 

-0.47732 

0.0 

mini 

-185.57294 

0.03494 

-0.53800 

-38.1 


a ICCI +Q energy in Eh- 


b zero-point energy in Eh, obtained from the harmonic CASSCF frequencies. 
c First column ICCI +Q energy plus zero-point energy + 185. Second column 
relative energy in kcal/mol. 
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Table IV. Computed stationary point geometries a b 



mini 

spl 

min2 

min3 

sp3 

min4 

sp4 

r 

1.029 

1.308 







1.003 

1.007 

1.045 

1.046 

0.991 

1.052 

1.124 

r N t Ni 

1.344 

1.260 

1.255 

1.253 

1.193 

1.247 

1.175 

*n 2 o 

1.201 

1.258 

1.404 

1.423 

1.606 

1.440 

1.941 

rw,o 


1.376 

0.983 

0.976 

0.981 

0.977 

0.986 

rHjO 







1.633 

Z H 1 N 1 N 2 

113.9 

81.2 






z H 2 N!N 2 

112.2 

118.5 

104.4 

103.3 

174.5 

108.7 

104.7 

Z N!N 2 0 

113.5 

104.5 

110.5 

107.8 

109.3 

110.5 

88.0 

Z NjHjO 


95.7 






Z HjONz 


78.6 

104.4 

101.5 

100.3 

102.4 

139.4 


fa c 

20.7 

fa d 

155.4 


° bond lengths in A and angles in degrees. 

b mini and spl are from an 8/6 CAS, while the remaining results are from a 10/10 

CAS. 

c angle between the H1N1N2 and N1N2O planes. 
d angle between the H2N1N2 and NiN 2 0 planes. 


17 



Table V. Computed frequencies and rotational constants (cm *) a . 



mini 

spl 

min2 

min3 

sp3 

min4 

sp4 

Wl 

3871 

3791 

3660 

3757 

4030 

3753 

3601 

u; 2 

3490 

2160 

3313 

3296 

3688 

3216 

2100 


1792 

1642 

1623 

1651 

1793 

1650 

1793 

U>4 

1709 

1562 

1466 

1462 

1194 

1454 

1174 

^5 

1368 

1348 

1405 

1366 

595 

1346 

801 

U; 6 

1207 

1006 

912 

884 

389 

838 

473 

Ct?7 

712 

1207 

652 

637 

576 

624 

998 

U?s 

647 

629 

959 

970 

344 

997 

330 

u>9 

542 

1990 i 

509 

364 

1359 i 

353 

1150 i 

A 

2.669 

2.416 

2.504 

2.737 

2.455 

2.318 

1.660 

B 

0.436 

0.504 

0.416 

0.400 

0.356 

0.415 

0.361 

c 

0.378 

0.417 

0.357 

0.349 

0.311 

0.352 

0.296 


“ mini and spl are from an 8/6 CAS, while the remaining results are from a 10/10 

CAS. 
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Table VI. Comparison of energetics with other calculations. 


present work 



ICCI+Q + ZPE 


BAC-MP4 C 

ICCI+Q^ 

HMW e 

CG-f 


refl a 

ref2 6 





NH 2 + NO 

0.0 

0.0 

0.0 

0.0 

0.0 


mini 

-44.0 

-38.1 

-48.1 

-45.7 

-42.3 

-28.9 

spl 

-14.4 

-8.5 

-19.7 

-13.3 

-5.6 


min2 

-45.7 

-39.8 

-47.5 

-47.9 

-42.1 

-31.3 

min3 

-44.5 

-38.6 

-45.5 

-46.5 

-40.3 

-30.4 

sp3 

-7.4 

-1.5 

-7.0 

-6.9 

2.8 


min4 

-45.5 

-39.6 

-47.2 

-47.3 

-41.7 

-31.7 

sp4 

-21.9 

-16.1 

-25.2 

-19.5 

-8.9 



a Energies computed with respect to H 2 + N 2 O and the NH 2 + NO asymptote 
positioned from experiment. 

b Energies computed with respect to NH 2 + NO. 
c Ref. 20. 

d Present results computed with respect to NH 2 + NO and without ZPE for com- 
parison to the work of HMW. 
e Ref. 22. 
f Ref. 18. 
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Figure Captions. 


Fig. 1. The electronic structure of mini, min2, min3, min4, and N 2 + H 2 O (prod). 

Fig. 2. Schematic representation of the minimum energy path for NH 2 + NO. The 
energies are from ICCI calculations (See the text.). 

Fig. 3. Geometries for the stationary points for reaction (1) and reaction (2). Mini 
is non-planar but the remaining structures are planar. 
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Theoretical Characterization of the Reaction 


NH2 + O — > products : 


Stephen P. Walch“ 
ELORET Institute 
Palo Alto, Ca. 94303 


Abstract. 

The potential energy surface for NH 2 + 0 has been characterized using complete 
active space SCF (CASSCF)/ derivative calculations to determine stationary point 
geometries and frequencies followed by internally contracted configuration interac- 
tion (ICCI) calculations to determine the energetics. The calculations predict an NO 
bond strength of 85.8 kcal/mol for NH 2 0. The barrier for isomerization of NH 2 0 
to trans-HNOH is predicted to be 48.0 kcal/mol and the barriers for H + HNO 
for mi ng NH 2 0 and NHOH are predicted to be 2.1 and 8.3 kcal/mol, respectively 
(all corrected for zero-point energy). 

The computed heats of formation for NH 2 0 and cis- and trans-HNOH reported 
by Soto, Page, and McKee are in good agreement with the present results. Our 
results also agree with those obtained by Melius and Binkley to within their stated 
± 3 kcal/mol error limits. 

The barrier for H + HNO — > H 2 4- NO is computed to be « 0.3 kcal/mol which 
is about 0.7 kcal/mol lower than the value computed by Soto and Page (without 
correction for zero-point energy). 

“Mailing Address: NASA Ames Research Center, Moffett Field, CA 94035. 
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I. Introduction 

NH 2 and NH play important roles in the oxidation of nitrogen containing species 
as well as in the reduction of NO e.g. in the thermal de-NO z process [1-4]. The 
reaction 


H + HNO -> H 2 +NO (1) 

may also be important in simulations of hypersonic flow, where heated air contains 
NO due to the thermal Zeldovich mechanism, and the reaction 

H + NO ^ HNO (2) 

followed by reaction (1) constitutes a pathway for recombination of H 2 with NO as 
a third body. 

Recently Dagdigian and coworkers [5-6] have studied the internal energy distri- 
bution in the products for the reactions, 


nh 2 + 0 - 

NH + 0H 

(3a) 

nh 2 + 0 - 

■* H + HNO 

( 36) 


which occur on the same potential energy surface (PES) as reaction (1). In both 
cases the OH and HNO products show less internal excitation than predicted by 
statistical theories. These authors suggested that additional ab initio calculations 
were needed to help interpret the experimental results. 

There have been a number of previous computational studies of the NH 2 O PES. 
Most have focused on the geometry and energetics of the stable isomers of NH 2 O. 
Melius and Binkley [7-8] used the bond additivity corrected M0ller-Plesset pertur- 
bation theory through fourth-order (BAC-MP4) method to study the NH 2 O and 
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HNOH (trans) isomers as well as some of the barriers to interconversion between 
the isomers of NH 2 O and the barrier to the abstraction reaction, which yields the 
same products as reaction (3a), but occurs on the quartet PES. More recently Soto, 
Page, and McKee [9] (SPM) have carried out a careful study of the 2 A r and 2 A" 
states of H 2 NO and the cis and trans isomers of HNOH, using Hartree-Fock (HF) 
plus configuration interaction (Cl) wavefunctions with basis sets through correla- 
tion consistent polarized valence triple-zeta (cc-pVTZ) [10]. Soto and Page [11] (SP) 
also studied the minimum energy path for reaction (1) using a cc-pVDZ (double- 
zeta) basis set [10] and carried out calculations of the rate constant as a function of 
temperature using canonical variational transition-state-theory. Since completion of 
this work, we also became aware of recent calculations by Page and Soto [12] which 
characterized the barriers for addition of H to HNO to give H 2 NO and HNOH. 

In the present work, these earlier studies are extended to include additional saddle 
points connecting the various NH 2 0 isomers and the possible product channels. The 
technical details of the calculations are given in Section II, the results are presented 
in Section III, and the conclusions are given in section IV. 

II. Computational Details. 

Two different basis sets were used in this work. For the complete active space SCF 
(CASSCF) derivative calculations the polarized valence double-zeta set (PVDZ) of 
Dunning and Hay [13] was used. The basis set for N and 0 is a (9s5p)/[3s2p] 
basis augmented by a single set of 3d functions with exponents of 0.80 and 0.85 
for N and 0, respectively. The H basis is (4s)/ [2s] augmented with a single set of 
2p functions with exponent 1.00. The basis set used in the Cl calculations is the 
Dunning correlation consistent triple-zeta natural orbital basis set (cc-pVTZ) [10] . 
This basis is [4s3p2dlf]. for N and O and [3s2pld] for H and is described in detail 
in Ref. 10. 
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The CASSCF/derivative calculations had 11 electrons distributed among 9 or- 
bitals (denoted by 11/9). Thus, all but the 0 Is, N Is, and O 2s like orbitals are 
included in the active space in the CASSCF calculations. In the subsequent inter- 
nally contracted configuration interaction (ICCI) calculations, the N Is and O Is 
electrons were not correlated. These Cl calculations were based on the 11/9 active 
space. The reference space for the ICCI calculations was restricted to no more than 
two electrons in the weakly occupied CASSCF orbitals. 

The CASSCF/derivative calculations used the SIRIUS/ABACUS system of pro- 
grams [14], while the ICCI calculations were carried out with MOLPRO [15-18]. 
A multi-reference analog of the Davidson’s correction [19] was added to the ICCI 
energies and is denoted by +Q. 

III. Discussion 

Table I gives the ICCI and ICCI+Q energies, the zero-point energies (based on 
the harmonic CASSCF frequencies), and the relative energies with respect to NH 2 
+ 0 and H + HNO. The energetics are also shown schematically in Fig. 1, while 
Fig. 2 shows the structures. The geometric parameters for each of the structures are 
given in Table II and the CASSCF harmonic frequencies and rotational constants 
are given in Table III. 

The NH 2 0 minimum energy structure in the present calculations has the NO 
tilted by 35.6 degrees with respect to the NH 2 plane. This is very similar to what 
was found by SPM, who obtained values of 41.1 and 26.0 degrees for the same angle 
for HF and Cl methods, respectively. As expected, the computed CASSCF bond 
lengths are found to be longer than the HF bond lengths reported by SPM. The 
computed L HNH is 114.0 0 in the present work compared to 115.6° and 117.5° 
obtained by SPM for HF and Cl, respectively. Calculations were also carried out 
for a planar structure, which is found to be a saddle point. From Table I it is 
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seen that the non-planar structure is slightly 0.6 kcal/mol) below the planar 
structure before correction for zero-point energy, but with inclusion of zero point 
energy the planar structure is 0.1 kcal/mol below the non-planar structure. Thus, 
in a one-dimensional model of the vibrational problem as a double-well minimum, 
even for the lowest vibrational level the zero-point energy is sufficient to surmount 
the barrier. For higher vibrational levels the vibrational energy is well above the 
barrier. This is consistent with the work of Pauzat et al. [20], where it is found 
that the vibrational levels for the NH 2 O out-of-plane mode look much the same for 
a flat potential as for a double-minimum potential with a small inversion barrier. 

The NH 2 O minimum is computed to be 57.1 kcal/mol below H + HNO. This 
compares to 58.3 kcal/mol obtained by SPM. Melius and Binkley [8] find NH 2 O to 
be 88.5 kcal/mol below NH 2 + 0. Using the computed separation between H + HNO 
and NH 2 0 from Table I plus the experimental separation between NH 2 + 0 and H 
+ HNO leads to 85.8 kcal/mol as our best estimate of the latter energy difference 
(vide infra). Since the BAC-MP4 calculations include some rather large empirical 
corrections, this level of agreement is probably satisfactory and the present results 
are expected to be more reliable. 

Table IV compares energetics obtained at the CASSCF level with those obtained 
at the ICCI level using several basis sets. As expected the CASSCF energetics do 
not agree well with the ICCI energetics. The significant point to be made with 
respect to Table IV is the comparison between the computed and experimental 
values of the locations of the NH + OH and H + HNO asymptotes. Here it is seen 
that the separation between NH 2 + 0 and NH + OH is accurately computed ( 
error of 1.5 kcal/mol) with the cc-pVTZ basis set, but even with a cc-pVQZ basis 
set the separation between NH 2 + O and H + HNO is in error by 5.3 kcal/mol. 
This is expected since the former separation involves breaking an NH bond and 
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forming an OH bond. The OH bond is slightly more difficult to describe in the 
calculation but the errors in the two bond strengths approximately cancel. The 
latter separation involves breaking an NH bond and forming an NO multiple bond. 
Clearly the NO bond is much more difficult to describe than the NH bond and the 
error in the separation reflects an expected error of about 2 kcal/mol in the NH 
bond but about 7 kcal/mol in the NO bond. Another way of saying this is that the 
former separation is isodesmic, while the latter is not. This argument justifies the 
procedure used above for computing the energy of NH 2 O with respect to NH 2 + 0 
as the computed separation between NH 2 O and H + HNO plus the experimental 
separation between H + HNO and NH 2 + 0. 

The structure spl is a saddle point interconnecting the NH 2 O minimum with 
trans-HNOH. The barrier is 48.0 kcal/mol with respect to NH 2 O. Thus, this saddle 
point is below the NH 2 + 0 energy. This barrier height compares well with the 
value of 50 kcal/mol obtained by Melius and Binkley [8]. 

From Table II it is seen that trans-HNNO is 5.2 kcal/mol below cis-HNNO. This 
is precisely the value obtained by SPM using the same basis set as used in the 
present calculations but with a Cl based only on the HF reference configuration. 
The computed binding energies with respect to H + HNO are 52.1 kcal/mol (trans) 
and 46.9 kcal/mol (cis). These compare very well with values of 52.2 kcal/mol 
(trans) and 47.0 kcal/mol (cis) obtained by SPM. 

Saddle points sp2 and sp3 correspond to adding H to the O end of HNO to give 
trans-HNOH and to the N end of HNO to give NH 2 O, respectively. The barrier 
heights with respect to H + HNO are 5.9 kcal/mol and 0.4 kcal/mol, respectively, 
for ICCI calculations at the CASSCF saddle points. Page and Soto [12] carried out 
very similar calculations, which gave 5.8 kcal/mol and 0.2 kcal/mol, respectively, for 
MRCI calculations using a cc-pVTZ basis set at the CASSCF saddle point. How- 
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ever, they found that reoptimizing the saddle point using MRCI increased the barri- 
ers by 2.4 and 1.7 kcal/mol respectively. This leads to 8.3 kcal/mol and 2.3 kcal/mol 
as our best estimates of the barrier heights for forming trans-HNOH and NH 2 O, 
respectively. Bozzelli and Dean [21] indicate a barrier of ~ 10 kcal/mol for sp2 
and a lower barrier for sp3 (taken from the work of Melius and Binkley [8]). They 
argue, based on analogy to other H addition reactions, that these barriers should 
be lower and use 3.5 kcal/mol for both sp2 and sp3. Our calculated results also 
indicate that the 10 kcal/mol barrier for sp2 obtained by Melius and Binkley is too 
large, but should not be as small as 3.5 kcal/mol. Also our calculations indicate 
that a barrier height of 3.5 kcal/mol is too large for sp3. Bozzeli and Dean used 
a barrier of 50.5 kcal/mol for spl compared to 48.0 kcal/mol in the present results 
and barriers of 3.5 kcal/mol for sp2 and sp3 and obtained a branching ratio of 14 % 
NH + OH (3a) and 86 % H + HNO (3b) in good agreement with the experimental 
result by Dransfield et al. [22]. 

SPM speculate that there could be a saddle point connecting NH 2 O or HNOH to 
H 2 + NO with energy below NH 2 + 0. The structure denoted by sp5 was obtained 
in an attempt to locate a saddle point for the process: 


cis - HNOH -> H 2 + NO (5) 

From Table III it is seen that sp5 has two imaginary frequencies, thus it is not a 
saddle point. The larger one corresponds to dissociation to H 2 + NO, while the 
smaller one corresponds to removing the plane of symmetry. (Note that this calcu- 
lation was carried out in Cj symmetry.) This structure is found to be 8.0 kcal/mol 
above NH 2 + O. This suggests that there is not a low energy pathway for reaction 
(5). 
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An attempt was also made to locate a saddle point for 


NH 2 0 -> H 2 + N0 (6) 

The inital geometry for the search was obtained by rotating the H 2 moiety of sp5 
about the H 2 center of mass to an orientation perpendicular to the plane. This 
search lead to sp3. The failure to find a saddle point for this process, as well as 
simple qualitative considerations of the orbital changes in this reaction, suggest that 
there is not a low energy pathway for reaction (6). 

The sp4 saddle point corresponds to reaction (1). The saddle point geometry 
which was obtained here using an 11/9 CASSCF has a longer th 1 H 2 than that 
obtained by SP using a smaller 5/5 CASSCF ( r of 1.41 A in the present 
work compared to 1.21 A obtained by SP). SP observed an elongation of th^h^ 
for Cl as compared to CASSCF. The same effect is seen in our results. Table V 
shows CASSCF and ICCI results (without inclusion of zero-point effects) along the 
minimum energy path (MEP) from sp4 to H + HNO. This MEP was obtained 
by following the CASSCF gradient away from a starting geometry where the sp4 
geometry was displaced slightly toward H + HNO and an ICCI calculation was 
carried out at each step in the walk. These results axe also shown in Fig. 3. From 
Table V it is seen that the barrier height on the ICCI surface is ~ 0.3 kcal/mol 
and corresponds to an of ~ 1.6 A. This barrier is significantly lower than 

the barrier of « 1.0 kcal/mol obtained by SPM using MRCI (without inclusion of 
zero-point effects). The difference presumably arises from the larger basis set used 
in the present work. 

SP predicted that the rate constant for reaction (1) should be about an order 
of magnitude larger than current experimental estimates. The lower barrier height 
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obtained in the present calculations would presumably mean that the computed 
rate is even larger, further increasing the discrepancy with experiment. 

IV. Conclusions. 

Calculations have been presented for the NH 2 + O potential energy surface us- 
ing complete active space SCF(CASSCF)/ derivative methods to characterize the 
stationary points and internally contracted configuration interaction (ICCI) to de- 
termine the energetics. 

The results for the 2 A' state of NH 2 O and for cis- and trans-HNOH are in good 
agreement with the work of Soto, Page, and McKee. In particular our results agree 
to within ± 1 kcal/mol with the heats of formation which they reported. 

The agreement with the BAC-MP4 results of Melius and Binkley is poorer. They 
appear to overestimate the binding energy in NH 2 O by ~ 3 kcal/mol and overesti- 
mate barrier heights by ss 2 kcal/mol. However, these errors are within their stated 
error bars of ± 3 kcal/mol. 

We predict barrier heights (zero-point corrected) of 2.1 and 8.3 kcal/mol for H 
-f HNO going to NH 2 O and NHOH, respectively. We also predict a barrier height 
(zero-point corrected) of 48.0 kcal/mol for isomerization of NH 2 O to trans-HNOH. 
This barrier height is less than the NH 2 -O bond strength which is computed to be 
85.8 kcal/mol. 

Production of H 2 + NO from NH 2 0 or cis-HNOH is predicted to involve barriers 
which are above the NH 2 4- O asymptotic energy. 

The barrier for H + HNO — ► H 2 + NO is predicted to be « 0.3 kcal/mol (before 

0-7 

zero-point correction). This is « T kcal/mol lower than the barrier used by Soto 
and Page in a canonical variational transition theory estimate of the rate of this 
reaction. Thus, the true rate may be even larger than predicted by Soto and Page, 
which in turn is about an order of magnitude larger than the experimental rate. 


9 



Thus, we concur with Soto and Page that a new experimental determination of this 
rate is needed. 


Acknowledgement. SPW was supported by NASA cooperative agreement num- 
ber NCC2-478. 
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Table I. Computed energies and zero-point corrections. 



Energy 0 

ZPE 6 

A E° 

A Ej 

nh 2 + 0 

-130. 74229(-. 76506) 

0.01848 

0.0 

28. 7 e 

NH 2 O(minimum) 

-130. 87396(-. 90028) 

0.02659 

-79.8 

-57.1 

NH 2 0(planar) 

-130. 86878(-. 89928) 

0.02545 

-79.9 

-57.2 

spl 

-130.78996(-. 81750) 

0.02021 

-31.8 

-9.1 

NHOH(trans) 

-130. 86597(-. 89218) 

0.02645 

-74.8 

-52.1 

NHOH(cis) 

-130. 85694(-. 88319) 

0.02569 

-69.6 

-46.9 

sp2 

-130. 76365(-. 78955) 

0.01612 

-16.8 

5.9 

sp3 

-130. 77316(-. 79753) 

0.01547 

-22.3 

0.4 

NH + OH 

-130. 75192(-. 77598) 

0.01539 

-8.8 

13.9 

sp4 

-130.77216(-. 79622) 

0.01184 

-23.7 

-1.0 

sp5 

-130.72053(-. 74861) 

0.01472 

8.0 

30.7 

H + HNO 

-130.77387(-. 79600) 

0.01330 

-22.7 

0.0 


a Energy in Eh obtained with the cc-pVTZ basis set. The first energy is the 
ICCI energy, while the energy in- parenthesis includes a multi-reference Davidson 
correction and is with respect to -130. Eh- 
b zero-point energy in Eh* 

c relative energy in kcal/mol with respect to NH 2 + O. 
d relative energy in kcal/mol with respect to HNO + H. 
c from experiment. 
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Table II. Computed CASSCF VDZP geometries®. 



nh 2 o 

spl 

NHOH(t) 

NHOH(c) 

sp2 

sp3 

sp4 

sp5 

*NH 1 

1.037 

1.049 

1.049 

1.053 

1.068 

1.079 

1.150 

1.289 

*nh 2 

1.037 

1.204 




1.924 

2.564 

1.264 

r no 

1.298 

1.448 

1.396 

1.393 

1.265 

1.230 

1.211 

1.351 

r oh 2 


1.214 

0.979 

0.980 

1.537 



1.408 

t H 1 H 2 







1.414 

0.963 

Z H!NH 2 

114.0 

121.8 







Z H 2 ON 


52.9 

102.2 

108.5 

113.3 




Z ONH 2 


53.7 







z NH 2 0 


73.5 







Z HiNO 


104.4 

99.9 

105.0 

105.9 

109.5 

109.1 

109.3 

Z H 2 NO 






119.6 



Z HaHjN 







174.5 


Z NHiH 2 








66.5 

Z H!H 2 0 








129.7 

Z H 2 ON 








54.5 

a 6 

35.6 








r 


138.6 



131.9 

116.2 




“ Bond lengths in A and angles in degrees. 
b angle between the NO bond and the bisector of Z H 1 NH 2 . 
c angle between the Hi NO and H 2 NO planes. 
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Table III. Computed CASSCF VDZP frequencies and rotational constants (cm J ). 



nh 2 o 

spl 

NHOH(t) 

NHOH(c) 

sp2 

sp3 

sp4 

sp5 

U?1 

3291 

3242 

3721 

3675 

2961 

2800 

1726 

2380 

U > 2 

1656 

2545 

3261 

3191 

1493 

1534 

1524 

1864 


1307 

1305 

1568 

1505 

1338 

1504 

1198 

1297 

u; 4 

691 

975 

1266 

1321 

697 

573 

490 

919 

w 5 

3421 

803 

1066 

1058 

586 

379 

258 

2278 i 

a? 6 

1304 

2201 i 

731 

529 

2013 i 

1204 i 

1724 i 

750 i 

A 

9.94 

8.93 

9.50 

9.41 

6.48 

5.01 

3.46 

6.87 

B 

1.01 

0.98 

1.03 

1.02 

1.10 

1.07 

1.07 

1.11 

C 

1.01 

0.93 

0.93 

0.92 

0.99 

0.94 

0.82 

0.96 
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Table IV. Computed relative energies (kcal/mol). 



CASSCF 

ICCI 





VDZP 

VDZP 

cc-VTZ 

cc-VQZ a 

Experiment 

nh 2 + 0 

0.0 

0.0 

0.0 

0.0 

0.0 

NH 2 O(minimum) 

-57.3 

-73.2 

-79.8 



NH 2 0(planar) 

-56.9 

-73.8 

-79.9 



spl 

-10.0 

-25.9 

-31.8 



NHOH(trans) 

-53.9 

-68.4 

-74.8 



NHOH(cis) 

-48.2 


-69.6 



sp2 

-5.8 


-16.8 



sp3 

-15.7 

-18.0 

-22.3 



sp4 

-18.4 


-23.7 



NH + OH 

-2.1 

-7.8 

-8.8 


-10.3 

H + HNO 

-21.5 

-19.1 

-22.7 

-23.4 

-28.7 

sp5 



8.0 




° This is a [5s4p3d2f/4s3p2d] basis set obtained from the cc-pVQZ basis set ( Ref. 
10.) by omitting the g function on N and O and the f function on H. 
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Table V. Computed energetics for reaction (l) a . 


iNHt 

t H 1 Hi 

1.150 

1.414 

1.106 

1.493 

1.097 

1.562 

1.091 

1.629 

1.088 

1.694 


CAS 

ICCI+Q 

4.41 

0.14 

4.12 

0.24 

3.63 

0.33 

3.15 

0.33 

2.69 

0.28 


a Bond lengths in A and energies in kcal/mol relative to H + HNO 
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Figure Captions. 


Fig. 1. Schematic representation of the potential energy surface for NH 2 + 0. The 
energies are from ICCI -f Q calculations. 

Fig. 2. Stationary point geometries for the NH 2 + 0 surface. See also Table II. 

Fig. 3. Potential for H + HNO — ► H 2 + NO. The minimum energy path is from 
C ASSCF / gradient calculations. The ICCI +Q energies were computed at the ge- 
ometries at each step on the walk. 
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Abstract. Complete active space SCF/ internally contracted configuration interac- 
tion (CASSCF/ICCI) calculations using large atomic natural orbital (ANO) basis 
sets are reported for CH 3 + O 2 . Two potential energy surfaces are found to be 
important in the CH 3 + 0 2 reaction. In C a symmetry, the lower 2 A" surface cor- 
relates with CH 3 + 0 2 ( 3 S“ ) and connects to a bound CH 3 00 species with no 
barrier, but leads only to CH 3 0 + O products. A higher surface of 2 A' symmetry 
correlates with CH 3 + 0 2 ( 1 A s ) and leads to CH 2 0 + OH with a computed barrier 
of 13.7 kcal/mol (with respect to CH 3 + 0 2 ( 3 £“ ))• Even in lower symmetry, two 
surfaces are involved leading to a more complex model for this reaction than had 
been previously considered. 

I. Introduction. 

The reaction 


ch 3 + o 2 - 

-+ CH z O + 0 

(la) 

CH Z +0 2 - 

♦ CH 2 0 + OH 

(lb) 


is believed to be an important chain branching reaction in hydrocarbon ignition 
processes [ 1 ]. The experimental situation has been discusssed by Zellner and Ewig 
[2]. The various determinations of the rate coefficient for reaction (la) are in rea- 
sonable agreement giving activation energies between 26.1 and 31.1 kcal/mol, which 
is consistent with the endothermicity of 28.7 kcal/mol [3]. However, there has been 
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considerable disagreement on the rate coefficient for reaction (lb). The determi- 
nations of vary by almost two orders of magnitude between 1500 and 2000K. 
A recent determination [4] of the rate of reaction (lb) gives an activation energy 
of only 9.1 kcal/mol. Zellner and Ewig have carried out RRKM calculations using 
a best estimate of the potential surface parameters for reaction (1) and assuming 
a model in which there is a common CH3O2 intermediate which can decompose 
to CH3O + 0 or rearrange to give CH 2 0 + OH. They conclude that channel lb 
will dominate over channel la at all temperatures below 2800K. There still remains 
some uncertainty in the experimental situation, since Baldwin and Golden [5] failed 
to observe reaction (lb) in an experiment at 1220K in contrast to the results of Ref. 
4. 

Given the importance of this reaction and the remaining uncertainty in the exper- 
imental situation, we have carried out ab-initio calculations of the potential energy 
surface (PES) to better define the energetics, particularly the barrier to formation 
of CH2O + OH. Section II describes some qualitative features of reaction (1). The 
technical details of the calculations are given in Section III. Section IV describes 
the results and the conclusions are given in Section V. 

II. Qualitative Features. 

Fig. 1 shows the electronic structure of the 2 A' and 2 A ;/ surfaces of CH3 + 
02- The presence of a mirror plane requires that one CH bond be in the C-0-0 
plane. Since the energy to rotate the CH 3 group is small, this constraint does not 
affect the energy significantly, but the use of a mirror plane reduces the size of 
the subsequent ICCI calculations considerably and this approximation was made 
here. As discussed in Section IV, the imposition of a plane of symmetry results in 
some stationary points in C 3 symmetry which have small negative eigenvalues of 
the Hessian matrix corresponding to rotation of the CH 3 group. However, the effect 
on the energetics is small. 

From Fig. 1 it is seen that the 2 A" surface correlates with the X 3 E“ ground 
state of O2. On this surface the distal O has an in plane doubly occupied 2p like 
orbital. Thus, 1,3-hydrogen migration to form CH200H ( and ultimately H2CO + 
OH) is not favorable. Therefore the product is expected to be CH 3 0 + O. Note 
that CH3O can form CH2O + H with a barrier of 23.3 kcal/mol [6]. 
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The 2 A' surface correlates with CH3 plus O2 in the 1 A ff state. On this surface the 
distal O has a singly occupied in plane O 2p like orbital and 1,3- hydrogen migration 
to give CH 2 0 + OH is favorable. 

III. Computational Details. 

Two different basis sets were used in this work. For the CASSCF derivative 
calculations the polarized double zeta set of Dunning and Hay [7] was used. The 
basis set for C and 0 is a (9s5p)/[3s2p] basis augmented by a single set of 3d 
functions with exponents of 0.75 and 0.85 for C and O, respectively. The H basis 
is (4s) / [2s] augmented with a single set of 2p functions with exponent 1.00. The 
basis set used in the Cl calculations is the Dunning correlation consistent triple zeta 
double polarization atomic natural orbital basis set [8]. This basis is [4s3p2dlf] for 
C and 0 and [3s2pld] for H and is described in detail in Ref. 8. 

The CASSCF/ derivative calculations had 7 electrons distributed among 7 orbitals 
(denoted 7/7). The active electrons correspond to the electrons in Fig. 1. which are 
in singly occupied orbitals. In the subsequent internally contracted configuration 
interaction (ICCI) calculations, the C Is and 0 Is electrons were not correlated. 
These Cl calculations were based on the 7/7 CAS calculation. The reference space 
for the ICCI calculations had the restriction that no more than two electrons were 
permitted in the weakly occupied CASSCF orbitals. 

The CASSCF/derivative calculations used the SIRIUS/ABACUS system of pro- 
grams [9], while the ICCI calculations were carried out with MOLPRO [10-13]. 
A multi-reference analog of the Davidson’s correction [14] was added to the ICCI 
energies and is denoted by +Q. 

IV. Discussion. 

The computed energetics for the 2 A' and 2 A" surfaces are given in Tables I and 
II, respectively. The zero-point energies are obtained from the CASSCF harmonic 
frequencies ( without scaling). The relative energies are computed with respect 
to CH3 + 0 2 ( 3 E7) on the 2 A" surface and include a multi- reference Davidson’s 
correction and zero-point energy. In each case the ICCI calculations are carried out 
at the geometries obtained from the CASSCF/derivative calculations. For the CH3 
+ 0 2 asymptote a supermolecule calculation was carried out at geometries obtained 
at the CASSCF level ( for the 3 E“ state in the case of 0 2 ). The energetics are also 
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shown graphically in Fig. 2. 

From Fig. 2 it is seen that, on the 2 A" surface, CH 3 adds to O 2 ( 3 £“) with no 
barrier to form CH 3 OO. The calculated barrier is actually slightly negative (See 
Table II.), this indicates that the true barrier is nearly zero and probably occurs for 
larger r C o than obtained at the CASSCF level. The CH 3 OO species is bound by 
27.6 kcal/mol, while the CH 3 0 + O product channel is computed to be endothermic 
by 26.2 kcal/mol compared to 28.7 kcal/mol from experiment [3]. From Fig. 2 it 
is also seen that CH 3 adds to 0 2 ( 1 A ? ) with a 3.7 kcal/mol barrier, leading to a 
2 A' excited state of CH 3 OO, which is 19.7 kcal/mol above the 2 A" ground state. 
There is a barrier of 13.7 kcal/mol, with respect to CH 3 + 0 2 ( 3 £ s )? to forming the 
CH 2 0 + OH product channel, which is computed to be exothermic by 54.1 kcal/mol 
compared to 53.2 kcal/mol from experiment [3]. No minimum corresponding to 
CH 2 OOH was found on this surface, rather the minimum energy path away from 
the saddle point proceeds directly to CH 2 0 + OH. 

Formation of the CH 2 0 + OH product channel from 0 2 ( 3 £~) requires coupling 
between the initially formed 2 A" CH 3 OO species and the exit channel portion of 
the 2 A' surface leading to CH 2 0 + OH. Thus, this problem involves two potential 
surfaces, which is more complicated than the model used by Zelner and Ewig [2]. 
The crossing region presumably involves geometries like those involved in the saddle 
point for 1 , 3 -hydrogen migration. As discussed in Section II, on the 2 A" surface 
moving the hydrogen toward the distal 0 is a repulsive interaction. Therefore, if C 5 
symmetry is retained, the 2 A" surface should rise more rapidly than the 2 A surface 
and the crossing would be expected to be before the saddle point. (This pathway 
is denoted by p 2 in Fig. 2 , while the pathway leading to CH 3 0 + 0 is denoted 
by pi.) If symmetry breaking is allowed (by rotation of the CH 3 group), the 2 A' 
minimum energy path would be displaced toward slightly higher energy, while the 
2 A" minimum energy path would be less repulsive. Therefore, the crossing would 
occur further along the reaction coordinate in the absense of symmetry. In this 
picture, the bottleneck would be the barrier to CH 2 0 + OH formation on the 2 A' 
surface, which is 13.7 kcal/mol. This is larger than the value of ~ 9 kcal/mol 
obtained in the recent experiments of Salto et al. [4], and which was also used 
by Zellner and Ewig [2] in their RRKM calculations. The ab initio calculations 
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reported here should be accurate to th ± 2 kcal/mol, based on the accuracy of the 
computed heats of reaction for (la) and (lb). This implies that the true barrier is 
larger than that used by Zellner and Ewig [2]. An additional complication is the 
need to treat the crossing of the two surfaces. Both of these factors would suggest 
that the true rate for formation of channel lb is smaller than estimated by Zellner 
and Ewig [2]. 

The CH3 + O2 energy given in Table I corresponds to the O2 (*A 5 ) state. The 
computed excitation energy to the O2 ( 1 A S ) state is 24.6 kcal/mol, which may be 
compared to the T 0 value [15] of 22.5 kcal/mol. Table I also gives the energy of the 
CH20 + OH asymptote with respect to CH 3 + 0 2 ( 3 E~ ). Here it is seen that the 
computed separation is within 1 kcal/mol of experiment. 

There is also a reported heat of formation of CH3O2 of 4.1 kcal/mol at 298K [16]. 
Our calculated value at OK is 7.2 kcal/mol, or about 5.1 kcal/mol at 298K [17], 
which is in reasonable agreement. This is especially true since the C, minimum 
in this case has a small imaginary frequency corresponding to symmetry breaking 
which would give a slightly lower value of the computed heat of formation if the 
geometry were fully optimized (vide infra). 

Tables III and IV give computed vibrational frequencies and rotational constants 
obtained for stationary points on the 2 A' and 2 A" surfaces. As mentioned in Section 
II forcing C 4 symmetry results in stationary points in C a symmetry which have small 
imaginary frequencies corresponding to symmetry breaking. This is seen in the case 
of the 2 A' surface both for the CH3-OO saddle point and for the CH3OO minimum 
( W12 = 82 i and 235 i ). However, the saddle point for CH 2 0 + OH production is 
a true saddle point in C a symmetry with an imaginary frequency of 2278 i. Similar 
effects are seen for the 2 A" surface in which there is a small imaginary frequency 
of 192 i for the CH3 00 minimum. An estimate of the effect of symmetry breaking 
involving rotation of the CH3 group can be obtained from the results of calculations 
on CH 2 0H at the SCF level. Here there are two C 3 symmetry structures with 
imaginary frequencies of 445 i and 338 i. These are 1.9 and 2.4 kcal/mol above the 
low symmetry true minimum. Thus, these effects introduce uncertainties as large 
as 2 kcal/mol in the CH3OO minima and one of the CH3-OO saddle points. This 
would make these stationary points slightly deeper than calculated here, but would 
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have no significant effects on the rates for la or lb, which depend on the barrier 
heights of the saddle points, which do not exhibit these symmetry breaking effects. 
Another minor problem with the frequency analysis is that the CH 3 OO saddle point 
on the 2 A" surface has a nearly zero frequency. This corresponds to a very floppy 
torsional mode. 

V. Conclusions. 

Complete active space SCF/ internally contracted configuration interaction 
(CASSCF/ICCI) calculations using large atomic natural orbital (ANO) basis sets 
have been carried out for CH 3 + O 2 . The calculations were carried out in C 3 sym- 
metry, which requires that one CH bond be in the plane defined by the C atom and 
two 0 atoms. Some of the stationary points obtained in C, symmetry have small 
imaginary frequencies corresponding to symmetry breaking. However, by compar- 
ison to calculations on CH 2 OH these symmetry breaking effects would have only 
small (<l -2 kcal/mol) effects on the energetics and do not effect the saddle point 
leading to the CH 2 0 + OH product channel which is the main concern in this paper. 

Two potential energy surfaces are found to be important in the CH 3 + O 2 reac- 
tion. In C, symmetry, the lower 2 A” surface correlates with CH 3 -+- O 2 ( 3 E 9 ) an d 
connects to a bound CH 3 OO species with no barrier, but leads only to CH 3 O + 0 
products. A higher surface of 2 A' symmetry correlates with CH 3 + 0 2 ( J A S ) and 
leads to CH 2 O + OH with a computed barrier of 13.7 kcal/mol (with respect to 
CH 3 + 0 2 ( 3 E~)). Production of CH 2 0 + OH requires a surface hopping from the 
2 A" to the 2 A' surface. This crossing is predicted to occur before the saddle point. 
The saddle point for formation of CH 2 O + OH is predicted to be the bottleneck to 
reaction. The computed barrier is larger than the barrier height of « 9 kcal/mol 
used by Zellner and Ewig in their RRKM calculations for this system. Both the 
higher barrier found here and the need for surface hopping would lead to a lower 
probability for formation of CH 2 0 + OH than predicted by Zellner and Ewig. 

The computed heat of formation of CH 3 O 2 is 7.2 kcal/mol at OK and is esti- 
mated to be 5.1 kcal/mol at 298K in good agreement with an experimental value 
of 4.1 kcal/mol. 
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Table I. Computed energies and zero-point corrections ( 2 A ; surface). 



Energy® 

zero-point energy 6 

A E c 

Exp. 

ch 3 + 0 2 

-189. 77857(-. 83292) 

0.03314 d 

24.6 

22.5 

ch 3 -o 2 

-189. 77189(-. 83048) 

0.03678 

28.4 


ch 3 o 2 

-189.83686(-. 89447) 

0.04295 

-7.9 


ch 3 oo -> ch 2 ooh 

-1 89. 79296(-. 85440) 

0.03735 

13.7 


CH 2 0 + OH 

-189. 90756(-. 96144) 

0.03632 

-54.1 

-53.2 


a The energies are in Eh a nd are i n the form ICCI(ICCI+Q +189.). 
b Zero-point energy in Eh based on the CASSCF harmonic frequencies. 
c Energy in kcal/mol relative to CH 3 + O 2 ( 3 S“) based on the ICCI+Q energies 
plus zero-point energy. 

d The O 2 zero point energy is obtained as 1/2 w 3 E“ (CAS) - 1/2 (a; 3 E £f (Exp.) - 
lo 1 A g (Exp.)) 
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Table Ha. Computed energies and zero-point corrections ( 2 A" surface). 



Energy® 

zero-point energy 6 

A E c 

ch 3 + o 2 

-189. 78797(-. 87226) 

0.03330 

0.0 

0 

a 

CO 

1 

0 

10 

-189. 79599(-. 87394) 

0.03421 

-0.5 

ch 3 o 2 

-189.86187(-. 92649) 

0.04354 

-27.6 


a The energies are in E// and are in the form ICCI(ICCI+Q +189.). 

^ Zero-point energy in E h based on the CASSCF harmonic frequencies. 
c Energy in kcal/mol relative to CH 3 + O 2 ( 3 E“) based on the ICCI+Q energies 
plus zero-point energy. 


9 


Table lib. Computed energies and zero-point corrections ( 4 A" surface). 



Energy 0 

zero- point energy 6 

A E c 

Exp 

CH 3 + 02 

-1 89. 79849(-. 87421) 

0.03330 

0.0 

0.0 

ch 3 o + 0 

-189. 78571 (-.83822) 

0.03905 

26.2 

28.7 


a The energies are in Eh and are in the form ICCI(ICCI+Q +189.). 

6 Zero-point energy in Eh based on the CASSCF harmonic frequencies. 
c Energy in kcal/mol relative to CH 3 + O 2 ( 3 E“) based on the ICCI+Q energies 
plus zero-point energy. 
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Table III. Computed stationary point frequencies and rotational constants(cm *). 
2 A' surface. 



CH 3 -00 

CH 3 OO 

CH 3 OO 


3343 

3363 

3263 

W 2 

3082 

3027 

1885 

U)$ 

1517 

1592 

1563 

Ct ? 4 

1054 

1515 

1156 

^5 

899 

1213 

1070 

U ) 6 

557 

988 

727 

U?7 

265 

776 

565 

U>8 

700i 

421 

2278i 


3441 

3330 

3390 

^10 

1466 

1531 

1210 

a> n 

522 

1199 

1135 

^12 

82i 

235z 

431 

A 

1.4049 

1.4483 

1.1196 

B 

0.2492 

0.3705 

0.4533 

C 

0.2213 

0.3129 

0.3459 


CH 2 OOH 
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Table IV. Computed stationary point frequencies and rotational constants(cm 1 ). 
2 A" surface. 



CHj-00 

CH 3 00 

Wl 

3350 

3261 

U>2 

3099 

3053 

W 3 

1504 

1597 

U 4 

1084 

1511 

^5 

495 

1219 

u; 6 

248 

1041 

(jjj 

138 

850 


156t 

534 

(Jj 9 

3444 

3320 

UlO 

1449 

1532 

Wll 

203 

1196 

<*>12 

ss 0 

192z 

A 

1.4988 

1.6532 

B 

0.1720 

0.3661 

c 

0.1594 

0.3180 
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Figure Captions. 


Fig. 1. The electronic structure of the 2 A' and 2 A" surfaces of CH3 + O2. 

Fig. 2 . Schematic representation of the lowest two potential energy surfaces of CH 3 
+ O2. The 2 A" surface is denoted by a solid curve, while the 2 A' surface is denoted 
by a dashed curve. For the 2 A" surface two product pathways are indicated. The 
pathway leading to CH 3 0 + 0 is denoted by pi, while the pathway denoted by 
p2 leads toward CH2O + 0 and is the region of crossing between the 2 A" and 2 A' 
surfaces. 
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Abstract 

Ab initio calculations employing large basis sets are performed to deter- 
mine an accurate potential energy surface for H— H 2 interactions for a broad range 
of separation distances. At large distances, the spherically-averaged potential de- 
termined from the calculated energies agrees well with the corresponding results 
determined from dispersion coefficients; the van der Waals well depth is predicted 
to be 75±3 fiEk- Large basis sets have also been applied to reexamine the accu- 
racy of theoretical repulsive potential energy surfaces (25-70 kcal/mol above the 
H— H 2 asymptote) at small interatomic separations; the Boothroyd, Keogh, Martin, 
and Peterson (BKMP) potential energy surface is found to agree with results of 
the present calculations to within the expected uncertainty (± 1 kcal/mol) of the 

a Mailing address: NASA Ames Research Center, Moffett Field, CA 94035-1000 
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fit. Multipolar expansions of the computed H-H 2 potential energy surface are re- 
ported for four internuclear separation distances (1.2, 1.401, 1.449, and 1.7 do) of 
the hydrogen molecule. The differential elastic scattering cross section calculated 
from the present results is compared with the measurements from a crossed beam 
experiment. 


TWork supported by Cooperative Agreement NCC 2-478 of the National Aeronau- 
tics and Space Administration (NASA) with ELORET Institute. 

I. Introduction 

The transport properties of H atoms in H 2 molecules are required for stud- 
ies of the preliminary heating of the fuel and subsequent combustion processes in a 
hydrogen-burning supersonic ramjet engine for the proposed National Aero-Space 
Plane (NASP). The accuracy of these transport properties is limited by the H- 
H 2 potential energy surface. While this potential energy surface is irrefutably the 
best characterized polyatomic potential energy surface available [1-9], the potential 
at large H—H 2 separation distances r, needed to determine diffusion and viscosity 
coefficients, is far less well studied than for small separations (r < 4 do). In addi- 
tion, the potential energy surface for shorter interatomic separation distances is of 
considerable interest for the theoretical determination of reactive scattering cross 
sections. Furthermore, there are discrepancies between calculated and measured 
results [10-13]. For example, the state resolved differential scattering cross sections 
for the D+H 2 — ► DH+H reaction shows a large discrepancy between the calculated 
and measured ratios of the backward-scattered to the forward-scattered intensi- 
ties [11-13]; an analysis [11] of the difference in the cross sections led to questions of 
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the accuracy of the high-energy potential energy surface (at about 1 eV above the 
H+H 2 asymptote), particularly for bent configurations. In this work we reinvesti- 
gate the H 3 potential energy surface using a much larger basis set than those used 
in previous work [1-9]. Our goals are to characterize the surface for longer H-H 2 
separation distances and to calibrate the existing potential energy surfaces in the 
repulsive region (25 to 70 kcal/mol above the H-H 2 asymptote), which has been 
suggested [ 10 - 11 ] as the important region for comparisons of theoretical results with 
measurements from reactive scattering experiments. 

II. Existing Potential Energy Surfaces 

The best available potential energy surface is that of Boothroyd, Keogh, 
Martin, and Peterson (BKMP) [9], which is thought to be accurate to within about 
1 kcal/mol. This potential is based on 772 ab initio points from several different 
sets of calculations. In order to include the different sets of computed points in 
their fitting procedure, BKMP developed a semiempirical basis set correction to 
adjust the calculated energies onto an equal footing. Hence, the BKMP potential 
energy surface supersedes the earlier Liu-Siegbahn-Truhlar-Horowitz (LSTH) [3] 
and double many-body expansion (DMBE) [7] surfaces as it incorporates the ab 
initio energies on which these potentials were based and includes several hundred 
additional points. It should be noted, however, that while only the DMBE potential 
behaves properly in the vicinity of the conical interaction, the BKMP, LSTH, and 
DMBE potentials yield si mil ar results for the scattering calculations [10]. 

While 772 points are used in the fit, none of the ab initio energies corre- 
spond to H-H 2 distances greater than 4 no- The range 4.0 < r < 6.0 ao defines 
an interpolation region between the ab initio data and an empirical potential en- 
ergy function that is constructed from the spherically- averaged potential Vo(r) of 
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Gengenbach et al. [14] and the anisotropic Vj model potential of Dalgaxno, Henry, 
and Roberts [15]. The spherical (Born-Mayer-spline-van der Waals (BMSV) type) 
potential of Gengenbach et al. [14] was derived from their measurements of absolute 
integral cross sections for D-H 2 scattering using the dispersion coefficients, C« ob- 
tained by Langhoff and Karplus [15] and C% from the results of Margenau [17]. The 
uncertainty in the measured cross sections (large error bars at low collision ener- 
gies [14]) and the difficulty in obtaining a unique potential fit [18,19] (e.g., changes in 
the van der Waals attraction can be compensated for by shifting the repulsive wall) 
limi t the accuracy of their potential at large r. Gengenbach et al. [14] estimated 
that the uncertainty in the well depth of the BMSV potential is ± 30%. The uncer- 
tainty in the van der Waals region is illustrated by the difference between the BMSV 
potential and the Born-Mayer dispersion (BMD) and modified BMD (MBMD) po- 
tentials. The latter two potentials were derived by Torello and Dondi [20] from their 
measurements of the differential cross sections for D-H* scattering using a crossed 
beam experiment and the combined data from their measurement and that of Ref. 
14, respectively, and the values of the dispersion coefficients calculated by Meyer, 
reported in Ref. [20]. The potential fit to the differential cross section is likewise 
not unique [18], but, nevertheless, is expected to be more sensitive to the van der 
Waals region since the rainbow structure in the scattering cross section has been 
resolved by their measurements. 

We compare the available spherically-averaged potentials in Figs. 1 and 2. 
Clearly, there is significant variation; this is disconcerting because this uncertainty 
would not allow a definitive determination of the transport properties that are 
described above. The BKMP potential [9] closely follows the BMSV potential as 
expected, but is unrealistic at large r because it falls below the interaction energies 
determined from the dispersion coefficients. The DMBE potential is much too 
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repulsive because it was constrained in the van der Waals region by six, scaled ab 
initio energies [4] that are not accurate. In addition to questions as to the accuracy 
of the existing van der Waals spherically-averaged potential, there has been no 
accurate determination of the leading anisotropic term V 2 (which is important for 
determining the inelastic contribution to the scattering from rotational transitions 
at low collision energies or large impact parameters [21]), but the general functional 
form has been predicted (see Sec. V below). 

III. Methods 

The coordinate system used in this work (pictured in Fig. 3) is that em- 
ployed in our previous study [8]. It is specified by the internuclear separation 
distance for the hydrogen molecule , the distance from the H atom to the center 
of mass of the H 2 molecule r and x ( tiie angle between a line from H to the H 2 
center of mass and the H 2 internuclear axis). Note that although refers to a 
molecule it is longer than t for some of our calculations. For convenience, selected 
calibration points were also performed using the coordinate representation used by 
BKMP [9], that is the coordinates r a , r 2 and 0 12 , where r 2 is H a -H* ? r 2 is H&-H c , 
and 0\2 is the angle between them. 

The determination of accurate potential energy surfaces for weakly interact- 
ing systems is in general a very demanding computational problem. The theoretical 
description must be able to describe the dispersion interaction and the basis set 
superposition error (BSSE) must be sufficiently small that it does not significantly 
affect the results. This puts severe demands on both the one- and n-particle basis 
sets. In this work we employ the multireference configuration interaction (MRCI) 
procedure used in our previous study [8] of H 3 that was shown to be essentially 
equivalent to a full configuration-interaction (FCI)j test calculations have demon- 
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strated that the MRCI energies are within 1 / lEk of the FCI energy. First a MRCI 
calculation is performed (in the self-consistent-field (SCF) molecular orbital basis), 
which includes all single and double excitations away from all distributions of the 
three electrons among 12 a' and 3 a" orbitals. This is followed by the same MRCI 
treatment using the natural orbitals (NOs) of the first MRCI. For all geometries 
the reference space comprises more than 99.9% of the final Cl wave function. 

While we use the same n-particle treatment, the previous basis sets [8] 
are unsuitable for the present study. The (5 s 4 p Id] basis set, which is essentially 
equivalent to the basis sets used to generate all of the ab initio points employed in the 
BKMP fit, is not flexible enough to describe the dispersion terms; it significantly 
underestimates the van der Waals well depth and is too repulsive in the region 
r = 3 to 6 a 0 . Even the larger [6s 5p 2d If] set is unsuitable because the H(H 2 ) 
BSSE is unacceptably large. Note that while these limitations make the basis sets 
unsuitable for the present study, they do not affect the previous results [8] to within 
their stated accuracy. 

Most of the current calculations are performed with a hydrogen one-particle 
basis set that is derived from the 10s set optimized for the 2 S state [22] supplemented 
with a diffuse s (0.02486). The innermost 6s functions are contracted based on the 
SCF orbital and the outer five functions are uncontracted. This is supplemented 
with the three p, two d, and one / function optimized by Dunning [23] and aug- 
mented with two diffuse p (0.1168 and 0.0467) and a diffuse d (0.2648) function, 
thus yielding a final basis set of the form (11s 5p 3d l/)/[6s 5p 3d 1/]. The sup- 
plementary diffuse functions are required to accurately describe the van der Waals 
region. However, the diffuse functions introduce linear dependency problems at the 
shorter r values. This occurs even for H 2 by itself at the shortest r values considered 
(0.8 and 0.9 ao). Thus, the linear dependency problems restrict, to some extent, 
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the geometries considered. For selected calibration calculations a diffuse / (0.32) is 
also added. The final basis sets employed replace the (3p 2d If) set with Dunning’s 
( 4 p 3 d 2/) H polarization set. This basis set is again augmented with two diffuse p 
(0.0984 and 0.0394) and a diffuse d (0.1972) function, giving a basis set of the form 
(11s 6 p 4 d 2/)/[6s 6 p Ad 2/]. The final basis set includes a diffuse / (0.35) and 
Dunning’s optimized g (2.358) giving a [ 6 s 6 p Ad 3/ 1 g\ basis set. Only the pure 
spherical harmonic components of the basis function are used. 

The basis sets employed in this work yield an energy for H that is only 
0.64 ixE h above the exact result. The H 2 potential energy curve (given in Table I) 
is in good agreement with the known non-relativistic limit [24] and the error in the 
H 2 dissociation energy is only 0.3 kcal/mol. The errors in the total energies are 
smaller than for the Slater basis used by Liu [l] to calculate the linear H 3 potential. 
The basis sets are sufficiently complete as to provide an accurate description of 
the dispersion terms. In addition, at t — 6.0 uoj the BSSE for all of the basis sets 
employed is less than 1 fiEh,- This is sufficiently small that the BSSE [H(H 2 )tH 2 (H)] 
is insignificant even in the van der Waals region; hence, we do not correct for BSSE 
in this work unless otherwise noted. 

IV. Calibration Points 

In the following discussion, we compare our computed results with the 
BKMP potential [9] and the data used for its construction; one major goal being to 
test the BKMP semiempirical basis set correction which allowed them to merge the 
ab initio energies from various groups in a consistent manner and to approximately 
correct for the known error in the H 2 interaction energies. The correction, Sb a3 , 
is obtained by using the London equation to provide a procedure for combining 
the H 2 basis set correction for the 3 interatomic distances in H 3 . On the basis of 
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test calculations, they estimated that the extrapolation energies were accurate to 
about 0.15 mEk for the linear symmetric conformations and about 0.3 mEh for 
the other conformations. In Table II we compare the basis set corrected energies 
compiled by BKMP from different authors with the results of this work. The ba- 
sis set correction for our computed points is about 2.5 times smaller than for the 
points computed by BKMP and about 1.7 times smaller than for the linear points 
of Liu. The agreement for the linear points, and particularly the symmetric points, 
is extremely good. Our extrapolated energies are slightly (0.07 kcal/mol) below 
the Monte Carlo energies [25], suggesting that the basis set correction is slightly 
overestimated for the linear points, but the difference is within the estimated uncer- 
tainty of the correction. For the nonlinear geometries, the agreement between the 
extrapolated energies decreases with increasing $12 * At 012—120°, our calculated 
energies are lower than BKMPs extrapolated values. The increase in the error in 
the BKMP extrapolated values with 0 12 is due to using a basis set correction that is 
independent of $12 for the case \ intuitively, however, the basis set correction 

must increase with the angle. 

Defining the correlation energy as the difference between the reference en- 
ergy and the MRCI energy (both in the natural orbital basis), we find that there is 
| more correlation energy at 120° than at 0°. Because of the very large basis sets 
used, it is logical to assume that the error in the correlation treatment is propor- 
tional to the correlation energy; this suggests that the basis set correction at 120° 
should be 33% larger than for 0°. We expect the underestimation of the correction 
to be true for all non-linear geometries; therefore the energies calculated in this 
work and those used to generate the BKMP potential are expected to differ the 
most for the large angle structures. 

In Table III, we compare a subset of the calculated energies that are more 
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than 25 kcal/mol above the IL-H 2 asymptote with the LSTH, DMBE, and BKMP 
potential energy surfaces. The root mean square (rms) differences for the 200 struc- 
tures that are more than 25 kcal/mol above the H-H 2 asymptote are 0.64, 0.49, and 
0.27 kcal/mol, respectively, for three potentials with the maximum differences being 
1.39, 1.94, and 0.85 kcal/mol. The London basis set correction relative to the H-H 2 
asymptote is also reported. As for the calibration points, the extrapolated energies 
are in excellent agreement with the BKMP potential for small angles. The basis 
set corrected values give an rms difference of 0.65, 0.54, and 0.34 kcal/mol with the 
maximum errors being 1.39, 2.15 and 1.02 kcal/mol. The largest errors occur for 
the large angle structures as previously discussed. The rms error for the 484 points 
computed in this work are 0.44(0.47), 0.37(0.40), and 0.19(0.23) kcal/mol for the 
three potentials with the basis set corrected results in parentheses. All the com- 
puted energies and the comparisons with the three potential functions are available 
as a PAPS document [26]. Selected cuts for a small value of r are compared with the 
corresponding results for the BKMP potential in Fig. 4; note that the agreement 
is very good and that the largest error lies within the claimed uncertainty of the 
BKMP potential. 

As disscussed by Liu [ 1 ] it is difficult to assess the accuracy of the calculated 
H 3 potential surfaces. The best procedure is to compare the calculated energy 
differences with improved treatments. As the basis set errors in the current work 
(either or the error in the H 2 energies) are about 40% smaller than those for the 
linear potential computed by Liu, we can assess the reliabilty of the data upon which 
the previous fits are based. Unfortunately, the same geometries were computed for 
only a few points. Following Liu’s [ 1 ] arguments, we expect our calculated H 3 
surface to lie between 0.1 and 0.5 kcal/mol above the true surface in the region 
1*0 < r < 4.0 a 0 . Relative to the calculated saddle point energies, the current 
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linear points appear to be within 0.1 kcal/mol of the energy separations computed 
by Liu over the specified range; the largest errors occuring for the larger r values. 
Comparing with the calculations of Siegbahn and Liu [2], the energy separations 
from the saddle point energy differ by as much as 0.5 kcal/mol for the large angle 
structures. While the errors are larger than the 0.1 kcal/mol estimated by Siegbahn 
and Liu [2], the errors are still relatively small. Also note, that the differences 
between our computed points and the BKMP fit (for the large angle structures) is 
about twice that compared with the differences in the ab initio calculations. We 
estimate that the current ab initio energies with respect to H+H+H are accurate 
to within about 0.5 kcal/mol, 1.5 times the error in our H 2 dissociation energy. 
The relative energies with respect to the saddle point should be considerably more 
accurate, probably to within 0.2 kcal/mol. Basis set incompleteness is still the 
principal source of error in the calculations; corrections for relativistic effects and 
Born-Oppenheimer corrections will be s m aller. 

V. Classical Barrier Height 

The classical barrier height is the energy difference between the saddle point 
and the H-H 2 asymptote. This quantity is an important aspect of the surface and 
has been extensively studied [1,5,8,25]. We calculate the barrier height to be 9.654, 
9.649, and 9.632 kcal/mol in the [6s 5p 3d 1/], [6s 6p 4 d 2/], and [6s 6p 4 d 3/ lg] 
basis sets, respectively (the [6s 6p 4 d 3/ lp] result is based on a single point cal- 
culation at n = r 2 = 1-75 a 0 ). Corrected for BSSE, the barrier heights are 9.670, 
9.663, and 9.643 kcal/mol. The computed barrier height may be compared with the 
9.65 kcal/mol calculated by Liu [5] (extended basis II) and is in very good agree- 
ment with the 9.59 ± 0.06 kcal/mol estimate of Liu [5]. The computed results are 
also in good agreement with the Monte Carlo results of Diedrich and Anderson [25] 
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of 9.61±0.01 kcal/mol where the error bar gives one standard deviation. The cal- 
culated total energies for the [6s 5 p 3 d If] and [6s 6p 4 d 2/] basis sets are 0.37 and 
0.26 kcal/mol respectively, above the Monte Carlo energies at r 2 = r 2 = 1.7 ao- We 
note that the reduction in the computed barrier height is not linear with the error 
in the H 2 energy so it difficult to extrapolate our results to predict the “true” bar- 
rier height. However, our results are consistent with Monte Carlo estimates. The 
London basis set correction for the [6s 5p 3d 1/] basis set is 0.12 kcal/mol giving an 
estimate of 9.53 kcal/mol; this is consistent with the observation that the London 
correction slightly overestimates the basis set correction. 

VI. H-H 2 potential energy 

The potential energy surface describing H-H 2 interactions can be conve- 
niently determined using the expansion 

= Vo (r, th, ) + V 2 (r, t H 7 )P 2 {cos x) + V 4 (r,r H2 )P 4 ( cos*) + ••• (1) 

where P 2n (cosx) is a Legendre polynomial. The ab initio points have been com- 
puted on a grid to facilitate the calculation of the expansion coefficients V 2n . The 
values of V 2 n were determined by a least squares fit of the calculated energies for 
fixed t and th 2 . We have tested the convergence of V 2 n and the stability of the fit 
using additional data points. The fitting error is small (< 0.02%) for r > 2 ao, but it 
increases rapidly at smaller r to a maximum value of 0.5% at r=1.6 uo- The results 
for V 2n are presented in Tables IVa-d. For consistency, the expansions (1) of the 
BKMP potential for comparison with the present results were generated with the 
same procedure as that described above. The leading anisotropic term of the ex- 
pansion (1) vanishes at the angle xo— cos 1 (l/''/3) =54.73561° for which P2( cos x) 
= 0; as would be expected from Tables IV, we find that the calculated values of 
the potential energy for this particular angle Xo yield a good approximation to the 
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spherically- averaged potential for a broad range of large r values. 

One expects that the potential energy deduced from the scattering mea- 
surements corresponds to an average separation distance for the hydrogen molecule 
that is about the same as the expectation value of the first vibrational state (i.e., 
t H2 = 1.449 a 0 ) rather than the equilibrium separation distance 1.401 a 0 . The 
highest collision energy of the scattering measurements by Gengenbach et al [14]. 
corresponds to about 34 kcal/mol; from Tables IV one finds that the spherically- 
averaged potentials for these two choices of differ by a few percent at the 
smallest value of r (~ 2 ao) probed by the experiment. On the other hand, the 
average potentials for these two choices are about the same for large r (> 4 a 0 ). 

Our calculated Vo is compared with other potentials in Figs. 1 and 2. The- 
agreement between the repulsive wall of our computed Vo and that of the BKMP 
potential is very good for small r. The value of 66 fJ*Ek for the van der Waals well 
depth D e of our calculated V 0 potential agrees very well with value 67 \iE k from the 
MBMD potential of Torello and Dondi [20]. Furthermore, our calculated Vo agrees 
well with the prediction of the asymptotic dispersion terms} for example, at t 
10 ao it is only 4 % lower than the slightly damped dispersion interaction energy (as 
described below). Thus overall, our computed spherically- averaged potential agrees 
very well with the best available data. We point out, in passing, that the repulsive 
wall of the LSTH potential in the van der Waals region agrees remarkably well with 
the results of the present calculation (see Fig. 2), notwithstanding the disclaimer 
stated in Ref. 3. 

Further improvement in the one-particle basis set will unquestionably yield 
a larger value for D e and also is expected to shorten r*, the position of the minimum 
of the potential well. It is difficult, a priori, to estimate the error in the calculations. 
The BSSE, which often can be used as a measure of basis set error, is negligible; 
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at r = 6.0 a 0 , the BSSE is less than 1 fiE h . Supplementing the basis set with a 
diffuse / function lowers Vo slightly, by 1 , 8 , 19, and 34 fiEh. at r=7.5, 6.0, 5.0, and 
4.0 <xo? respectively. Expanding the polarization set, going from the (5 p 3d If) to 
the ( 6 p 4 d 2 /) set has a slightly smaller effect; the ( 6 p 4 d 2 /) results are -0.1, 4.7, 
and 9.1 fiEh (at r=7.5, 6.0, and 4.0 a 0 ) above the diffuse / basis set results. The 
results indicate that the [6s 5 p 3d If] basis set employed in this work is reasonably 
complete even for the van der Waals region-the error in the H 2 energy will largely 
cancel. The largest deficiencies in the basis sets are the lack of diffuse / and g 
functions necessary to describe the higher-order polarizabilities. However, based on 
the calibration results we estimate that the true van der Waals well depth to be no 
more than about 20 % deeper than the present value 66 fiEh ; that is, we estimate 
D e to be 72 ± 6 fiE h - To improve upon this estimate, we have employed the larger 
[ 6 s 6 p Ad 3/ 1 g] basis set. Using the modifed potential energy fit derived below and 
our estimated D e of 72 fiEht we predict a bond length of 6.54 <io- A.t this predicted 
bond length, we obtain a D e value of 72.4 fiE h , which is unquestionably a lower 
bound (the BSSE is only 0.23 fiE h )- These results lead to our recommended value 
of 75 ± 3 fiEh.- From the foregoing discussion, we conclude that the values of both 
r e and D e from the BMSV potential are too large. 

As mentioned above, the anisotropic coefficient V 2 is of particular signifi- 
cance for rotational excitation by collisions. The general functional form was de- 
duced by Tang [27] based upon a semiempirical potential by Tang and Karplus [28] 
and the asymptotic long-range form. r H 2 ) f° r values of rn 2 near the equilib- 

rium distance has 3 sign changes; the long-range V 2 is negative and is determined 
by the anisotropy of the dispersion interaction. At large r, including the van der 
Waals region, as well as the region of the long-range dispersion forces, the inter- 
action potential is monotonically increasing from 0 to 90° (with r and th 2 fixed, 
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e.g., see Fig. 5). In the range 3.7 < r < 5.2 ao, Vz is positive indicating that a 
90° approach is the most favorable. At slightly smaller r values, Vz is negative 
and sizable (e.g., see Fig. 6) and a linear approach is strongly favored. This occurs 
because the repulsion between the H 2 and the H is smallest for this orientation; the 
quadrupole moment of H 2 is positive. At shorter range, r < 1.9 ao, Vz is positive 
again and increases sharply; this is to some extent an artifact of the coordinate 
system employed as the minimum H-H distance is now much smaller than th 2 for 
a linear approach. 

In Fig. 6, we compare V 2 and V 4 deduced from this work with the corre- 
sponding results from the BKMP potential. The agreement is excellent for small r; 
for example, the inner minimum of V 2 is only slightly more shallow than that derived 
from the BKMP fit. The ab initio values are more realistic (smoother), however, 
in the interpolation region at large r as would be expected. Improvements in the 
basis set are expected to have only a small effect on the anisotropic coefficients as 
the errors cancel in the H 2 energy differences. 

The ab initio results for the van der Waals region at large r are shown in 
Fig. 5. Note the uniform spacing of the potential energy curves with respect to x; 
this is consistent with equal differences from the leading term (note that AP 2 =3/8 
for adjacent angles) of the expansion (1) when the higher-order terms axe negligible. 

We have found that the potential energy in the van der Waals region can 
be described well by a modified version [29] of the potential form developed by Tang 
and Toennies [30] that can be written 

V(r,x) = V S r(t,x) + V D LR.(r,x) ( 2 ) 

where the short-range repulsive potential Vsr has the Born-Mayer form 

V S r{t,x) = A(x)exp[-a(x)r] (3) 
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and the long-range attractive dispersion energy is obtained from 

V DL R{r,x) = ( 4 ) 

71 “ 3 

using the universal form of the damping function that can be related to an incom- 
plete gamma function [31] 

2n x k 

f 2n (x) = 1 - exp(-a;)]r -£J. (5) 

Jk=0 

The parameters A and a -1 characterize the strength and range, respectively, of 
the repulsive interaction. The long-range forces are determined by the dispersion 
coefficients that are approximated by the expansion 

C 2n (x) = <?2n[l + r 2 niMcOSx)]. (6) 

The isotropic components C 2n calculated by Meyer (reported in Ref. [20]) are 8.813 
alE h , 163.87 a^Eh, and 4232.2 aj° E k for n = 3-5, respectively. The theoretical 
value of Ce agrees with the corresponding value deduced by Zeiss and Meath [32] 
from measured oscillator strengths, photoabsorption, and scattering data. Higher- 
order values can be obtained from the recursion relations [30,33] 

C 2n +4 = Qn(C 2n +2 / C 2n Y C 2 -n- 2 (7) 

with the values of Cl n calculated for hydrogen [33]. 

We have calculated the leading anisotropic contribution T 6 from the com- 
bining relations developed by Fuchs et al. [34] using polarizability data (exact results 
for H [35] and calculated/measured data for H 2 [36]) and the calculated dispersion 
coefficients for H-H and H 2 -H 2 interactions. The value of C 6 for H 2 has been de- 
termined exactly [33]; the values of C% aLbM for the molecular interactions were 
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obtained from a transformation [ 37 ] of the results for the angular representation 
£,L A L B L ca j cu ] a ^ e( ^ ijy Ptiyks and Wormer [38]. We obtain the result 1?6 = 0.100, 
which agrees well with the value 0.114 determined by Langhoff et al. [39] from op- 
tical dispersion, absorption, and scattering data. By comparing the contribution 
from the leading anisotropic dispersion term described above with the values of V 2 
determined from the ab initio data at large r, we can estimate the contribution from 
the next higher-order term. Thus we obtain the result Ts ~ 0.40 from the data for 
r in the range 8.5-10.0 ao. 

Following the discussion of Ref. 29, we can readily determine the repulsive 
parameters from the ab initio data for the parallel and perpendicular orientations 
(x = 0° and 90°, respectively); i.e., we obtain 

lnA(x) = 1-390 + 0.233P 2 (cos x) (8) 

and 

o:(x) = 1.629 + 0 . 018 P 2 (cosx). (9) 

We find that the result obtained with x = Xo from the modified Tang- 
Toennies (MTT) potential energy surface of this work (i.e., Eqs. (2-9) with the 
above data for the dispersion coefficients) agrees well with the spherically-averaged 
potential (i.e., V 0 of Table IVc). In addition, if we constrain the spherically- averaged 
potential to the ab initio results of Table IVc for the upper region of the repul- 
sive wall, the long-range behavior described by the dispersion interaction described 
above, and the expected value of D e = 75 fiEh from the above analysis; we obtain 
A = 5.05 fiEh. and a = 1.669 a q 1 . The value of r e from this predicted potential 
is 6.51 ao which is slightly less than the value of 6.54 ao deduced by Torello and 
Dondi [20] from their measured differential scattering data. The results for the pre- 
dicted upper and lower bounds for 1 D e are also reported in Table V with those for 
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the calculated and measured potentials described above. The MTT potential for x 
= xo Is compared with the calculated data and also the predicted potential with 
the deeper potential well in Fig. 7. 

As described above, the potentials obtained from the differential scattering 
measurements agree well with the present results in the region of the van der Waals 
minimum; the repulsive wall of these potentials is, however, considerably weaker 
(see Fig. 2). Because potentials fitted to scattering data may not be unique, it 
is of interest to compare the scattering from the present potential data with the 
measured results. 

We have calculated the differential elastic scattering cross section using the 
sudden approximations of Parker and Pack [40] from scattering phase shifts that 
have been determined by direct numerical integration of the Schrodinger equation. 
The potential energies for these calculations were constructed from spline fits to the 
ab initio data of Table III and the long-range tails from the MTT potential for the 
values of x specified in Fig. 5. In addition, we have also calculated the scattering 
from the spherically-averaged potential (i.e., the values of Vo from Table IVc and 
the long-range tails from the MTT potential with % = Xo)* The cross sections 
from these two different approaches are found to be nearly identical, as would be 
expected when the anisotropic terms of the potential energy (1) are small. 

The theoretical scattering results are compared with the experimental data 
of Ref. 20 in Fig. 8. We have averaged the cross section over the velocity profile 
of the primary beam, but the theoretical results represent a high-resolution cross 
section, in that, we have not taken the beam size, detector area, etc. into account; 
unfortunately, sufficient data was not provided by the authors to allow this averaging 
of the scattering. In light of this difference, the agreement between theory and 
experiment can be considered satisfactory. 
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VII. Conclusions 

Ab initio calculations have been performed using a very large basis set to 
calibrate the H-H 2 surface. The basis set incompleteness is much smaller than in 
all previous calculations. The results confirm the accuracy of the recent Boothroyd, 
Keogh, Martin, and Peterson (BKMP) potential and their scheme for correcting 
the H 3 calculations using the error in the H 2 potential. Our calculations do show, 
however, that the accuracy of the correction decreases with increasing angle, as 
some of our computed points are below their extrapolated energies. Our computed 
data fill in a region of the potential that has been left out of previous theoretical 
studies; that is, the region of large r from 4 ao out to about 8—10 do where the 
long-range dispersion interaction energies are sufficiently accurate. The differential 
elastic scattering cross section calculated from the present results for the potential 
energies agrees well with the measured data. We have derived a potential energy 
function for the van der Waals region at large r to augment the available potential 
energy fits that cover the small r region. The present results will allow an accurate 
determination of the contribution from H-H 2 interactions to transport properties; 
collision integrals for hydrogen will be reported in a future publication. 
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Table I. Total energies (Eh) for the X *11* state of H 2 . 


Basis set 


K a o) 

[6s 5 p 3 d 1/] 

[6s 6p 4 d 2 /] 

[6s 6p Ad Zf 

0.8 

-1.01964331“ 

-1.01964589“ 


0.9 

-1.08290665“ 



0.95 

-1.10569977 



1.0 

-1.12388052 

-1.12417145 


1.05 

-1.13825085 



1.1 

-1.14945650 



1.15 

-1.15802352 



1.2 

-1.16438212 

-1.16459398 


1.25 

-1.16888723 



1.3 

-1.17183355 



1.4 

-1.17399454 

-1.17416126 

-1.17424203 

1.401 

-1.17399505 



1.449 

-1.17360611 

-1.17376375 

-1.17384418 

1.5 

-1.17240057 



1.6 

-1.16815089 

-1.16828665 


1.7 

-1.16204450 



1.75 

-1.15848532 

-1.15860666 

-1.15868356 

1.8 

-1.15466972 

-1.15478735 


1.9 

-1.14646449 



2.0 

-1.13775776 

-1.13786438 


2.2 

-1.11977532 

-1.11987628 


2.4 

-1.10208144 

-1.10218103 


2.6 

-1.08546434 



2.8 

-1.07037029 



3.0 

-1.05702837 

-1.08556502 


3.4 

-1.03581724 



3.8 

-1.02134469 



4.0 

-1.01621290 



4.5 

-1.00787702 



5.0 

-1.00371123 



5.5 

-1.00172715 



6.0 

-1.00080524 

-1.00081273 


6.5 

-1.00038078 



7.0 

-1.00018488 



8.0 

-1.00004951 

-1.00005118 


9.0 

-1.00001649 



10. 

-1.00000661 
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12 . 

100 . 


- 1.00000106 

- 0.99999868 


- 0.99999868 


- 0.99999868 


a The Cl vectors in the MRCI procedure were picked based on overlap with the 
previous iteration to avoid “collapse” due to linear dependence in the basis set; the 
values might therefore be less accurate than the others. 
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Table II. Comparison of the basis set corrected II 3 energies. 
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Table III. Comparison of computed H3 energies with the LSTH, DMBE, and BKMP potentials. 
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the text. The basis set correction for II 2 is —0.3133 kcal/mol. 



Table IVa. Expansion coefficients' 1 V n (kcal/mol) for = 1.2 ao 


r(a 0 ) Vo V 2 V 4 V 6 V 8 


1.8 

48.477 

7.206 

2.0 

35.987 

1.861 

2.2 

27.216 

0.045 

2.5 

18.176 

-0.331 

3.0 

9.183 

0.059 

4.0 

1.934 

0.102 

5.0 

0.253 

0.009 

6.0 

-0.020 

-0.008 

7.0 

-0.033 

-0.006 

8.0 

-0.019 

-0.003 


5.113 

0.509 

0.102 

2.302 

0.237 

0.032 

0.993 

0.112 

0.012 

0.251 

0.035 

0.003 


0.016 

0.000 

- 0.001 

0.000 

0.000 

0.000 


“Blank entries indicate that the term was not included in the fit. 
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Table IYb. Expansion coefficients V^kcal/mol) for ry, = 1.401 ao 


’’( ao ) 

^0 

^2 

V 4 


V 8 

1.8 

43.386 

6.051 

10.473 

1.049 

0.430 

2.0 

31.794 

-1.320 

4.766 

0.530 

0.129 

2.2 

24.036 

-3.332 

1.995 

0.275 

0.043 

2.5 

16.322 

-2.793 

0.391 

0.093 

0.011 

3.0 

8.633 

-0.858 

-0.055 



4.0 

1.960 

0.048 

-0.010 



5.0 

0.266 

0.007 

-0.002 



6.0 

-0.023 

-0.011 

-0.001 



7.0 

-0.038 

-0.008 

0.000 



8.0 

-0.022 

-0.004 

0.000 
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Table IVc. Expansion coefficients V n (kcal/mol) for rjj 2 — 1.449 ao 


r ( a o ) 

Vo 

v 2 

V 4 

V 6 

V 8 

1.6 

60.2286 

26.9028 

25.7629 

2.9543 

1.9716 

1.7 

49.9297 

13.9850 

17.8003 

1.8156 

1.0662 

1.8 

41.9195 

5.7196 

12.2588 

1.2063 

0.5942 

1.9 

35.6059 

0.6223 

8.3590 

0.8430 

0.3225 

2.0 

30.5564 

- 2.3270 

5.6128 

0.6115 

0.1764 

2.2 

23.0717 

- 4.4369 

2.3519 

0.3294 

0.0570 

2.5 

15.7320 

- 3.6383 

0.4373 

0.1144 

0.0151 

3.0 

8.4351 

- 1.2047 

- 0.0880 

0.0092 

0.0027 

3.5 

4.2592 

- 0.1971 

- 0.0468 

- 0.0004 

0.0000 

4.0 

1.9579 

0.0239 

- 0.0145 

- 0.0002 

- 0.0004 

4.5 

0.7985 

0.0295 

- 0.0046 

- 0.0002 

- 0.0004 

5.0 

0.2684 

0.0061 

- 0.0020 

- 0.0002 

- 0.0003 

5.5 

0.0514 

- 0.0077 

- 0.0013 

- 0.0001 

- 0.0002 

6.0 

- 0.0240 

- 0.0115 

- 0.0010 

0.0001 

- 0.0001 

6.5 

- 0.0417 

- 0.0107 

- 0.0007 

0.0001 

0.0000 

7.0 

- 0.0389 

- 0.0083 

- 0.0005 

0.0002 

0.0000 

7.5 

- 0.0309 

- 0.0061 

- 0.0003 

0.0001 

0.0000 

8.0 

- 0.0230 

- 0.0043 

- 0.0002 

0.0001 

0.0000 

8.5 

- 0.0167 

- 0.0030 

- 0.0001 

0.0000 

0.0000 

9.0 

- 0.0120 

- 0.0020 

- 0.0001 

0.0000 

0.0000 

10.0 

- 0.0064 

- 0.0010 

0.0000 

0.0000 

0.0000 
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Table IYd. Expansion coefficients V n (kcal/mol) for th 2 — 1.7 ao 


Kao) 

^0 

v 2 

v 4 

^6 

v 8 

Vio 

1.8 

32.857 

4.962 

25.585 

2.265 

2.843 

0.270 

2.0 

22.657 

-8.380 

12.465 

1.083 

0.813 

0.029 

2.2 

16.657 

-11.898 

5.486 

0.665 

0.244 

0.013 

2.5 

11.563 

-9.961 

0.949 

0.308 

0.050 


3.0 

6.885 

-4.042 

-0.427 

0.019 

0.008 


4.0 

1.884 

-0.219 

-0.068 




5.0 

0.276 

-0.009 

-0.006 




6.0 

-0.029 

-0.017 

-0.002 




7.0 

-0.045 

-0.012 

-0.001 




8.0 

-0.027 

-0.006 

0.000 
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Table V. Summary of the results for the H-H 2 van der Waals potential well. 


r e (a„) D e {nE h ) A (E h ) a(a^) 


ab initio(rH 7 = 1.449 ao) 

6.62 

lower bound 

6.54 

predicted 

6.51 

upper bound 

6.48 

BMD a 

6.54 

MBMD a 

6.56 


67 

4.01 

1.629 

72 

4.55 

1.652 

75 

5.05 

1.669 

78 

5.57 

1.686 

61.4 



66.8 




“Ref. 20. 
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Figure Captions 


Fig. 1. A comparison of the spherically- averaged H-H 2 interaction energy for 
r H2 =1.401 a 0 from the present calculation (solid line) with the corresponding av- 
erages for various fitted potentials (dotted lines) and the potential energies deduced 
from measured scattering cross sections (dashed lines). 

Fig. 2. A comparison of the spherically- averaged H— H 2 interaction energy from the 
present calculation (solid line) with the corresponding averages for various fitted po- 
tentials (dotted lines) and the potential energies deduced from measured scattering 
cross sections (dashed lines). 

Fig. 3. The coordinate system used for most of the energy calculations. 

Fig. 4. Selected cuts through the potential energy from the present calculation 
(data points with a spline curve fit) obtained for r = 1.8 <Xo &ud certain values of 
r . The corresponding cuts for the BKMP potential (dashed line) are shown for 
comparison. 

Fig. 5. Van der Waals potential energy curves for r H2 = 1.449 a 0 and certain values 
of x\ the curves were constructed from spline fits to the discrete data. 

Fig. 6. A comparison of the potential coefficients from the present calculation (data 
points with a spline curve fit) for tjj 2 — 1.449 o-q with the corresponding results for 
the BKMP potential (dashed fines). 

Fig. 7. A comparison of the values of Vo(r,rj/ 2 = 1.449 ao) from Table IVc (data 
points) with the values obtained from the MTT potential for x = Xo (dashed line) 
and the predicted potential (solid fine) with D e = 75 which is like the MTT 

potential for x = Xo except that the repulsive parameters have the values indicated 
in the figure. 

Fig. 8. A comparison of the theoretical high-resolution D— H 2 differential elastic 
scattering cross section (solid fine) calculated from the present results for the in- 
teraction energy (see text) with the measured data (represented by the data points 
with error bars) from the crossed beam scattering experiment of Ref. 20. The cross 
sections are shown for the laboratory system for a peak collision energy = 35.8 meV. 
The measured cross sections have been normalized to the theoretical results at large 
scattering angles since absolute cross sections were not obtained in the experiment. 
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